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Benjamin Franklin performed a famous 


experiment 200 years ago, which had a 
counterpart in the investigations of ground 
potentials carried on much later by the 
Schlumbergers. Franklin proved the exist- 
ence of positive and negative electrical 
charges by flying a kite in a rainstorm 

This early experiment in the field of natural electrical potentials en- 
couraged other investigators to start and continue work of increasing prac 
tical value. The measurement of ground potentials in bore holes by the 
Schlumbergers was one such important development. The resulting sponta 
neous potential curve was added to the already existing resistivity curve to 
complete the fundamental Electrical Log. The development, now only twenty 
years old, has had an increasing influence on oil exploration and pro- 


duction 


neered the electrical 


log and has been 


in the field of electrical oil well 
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WOULD YOU RISK YOUR DRIVE TO SAVE A BALL? 


Sometimes you may get a pretty 
fair drive by teeing off with an old 
ball to save a new one. And sometimes 
you get a pretty fair oil well even 
though you attempt to perforate the 
casing with out-moded methods. But 
the record of thousands of JET per- 
forated wells shows that you best 
protect your well investment, and ob- 
tain increased productivity by per- 


forating with Welex Jets. That 


is why more than 50° of all cas- 
ing perforating done in the U. S. 
today is done by the JET process origi- 
nated by Welex. Welex Jet Perforating 
gives you the economy of positive 
penetration and known hole density 
... with no burr, no fracturing, and 
no debris left in the well. Results 
cost less with Welex Jets. Call your 
nearest Welex Station for prompt 


service day or night. 


Welex 


SET SERVICES INC. 


A GENERAL OFFICE: 3909 Hemphill Street © Fort Worth 9, Texas 
FIELD STATIONS; Shreveport @ Ardmore @ Lindsay @ Shawnee @ Corpus Christi 
Falfurrias @ Houston @ Kilgore @ Odessa © Hobbs @ Wichita Falls 





RESERVOIR ENERGY CAN BE SAVED 
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OIL PRODUCTION CAN BE INCREASED 





Series 
at the 





Any number of Drain Holes, devi- 
ating with a 16-18 ft. radius, at 
different levels and in different 
oriented directions from a well and 
EXTENDING for 60 feet or MORE 
horizontally into the oil sand, can 
be drilled from any old well or 
from any new well AT ANY DEPTH 
...- with Standard Rotary Drilling 
Methods and with conventional © 
4%" Rock Bits. 


In any type of Formation and Reservoir Conditions, oil 





will flow slowly from all sides into the Drain Holes, from 





which it flows slowly and without friction to the well. 


ZUBLIN DRAIN HOLES CONTROL THE FLOW OF OIL 
WITHIN THE OIL SAND SURROUNDING THE WELL 





Drain Holes are 
open and do not 
require casing for 
their protection. 
They are full of 
high pressure fluid 
at all times. 


PAINFUL PRODUCTION LOSSES of Oil and Gas 
Result from the following Causes: 


1. Tight Oil Sands with low permeability build up 
tremendous friction which hinders the flow of oil to 
the well. 


2. Loose Oil Sands with high permeability allow oil to 
flow at such high speed through the formation that 
the oil carries sands, colloidal and residual matters 
with it, depositing them either outside or inside the 
well with the following two-fold damages: 

Outside Deposits clog the sand pores of the formation 
and build up a filter cake around the well which 
gradually chokes the well to a fraction of its normal 
production. 

Inside Deposits sand up the well or its pumping equip- 
ment which then requires frequent pulling jobs and 
replacements with the corresponding loss of produc- 
tion, time and money. 

3. Improper Water Drive produces excess water due 
to coning around the well. Bottom water, with less 
resistance to flow than oil, fills the depleted sand 
around the well, which then produces water instead 
of oil. 

4. High Gas-Oil Ratios are created when the gas in 
solution rapidly becomes free gas due to agitation 
of the oil during its high speed—turbulent flow. This 
gas, by-passing the oil, reaches the well first, leaving 
the oil behind. 


NOW — please note that, in every instance, the 
loss of production is caused either by oil flowing 
too slow or oil flowing too fast. 


Thus, unbelievably large amounts of oil remain in 
the ground because there has been no planned, 
effective Flow Control of oil within the oil sand. 


ZUBLIN DRAIN HOLES have been invented to regulate 
the flow of oil and gas and to prevent or reduce the 
flow of sand and water within the formation BEFORE 
they reach the well. 














CALL JE 6151 or JE 4433 for a ‘Full Size’ Demon- 
stration in our plant (2369 East 51st St., Vernon, Los 
Angeles) showing how Drain Holes are drilled with 


ZUBLIN FLEXIBLE 
DRILL PIPE 


The Equipment is for rent to Producers and Contractors 
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SINCLAIR RESEARCH LABORATORIES—nine buildings containing the 
most modern testing equipment known—have contributed many of 
today’s most important developments in petroleum products, pro- 


duction and refining. Under the Sinclair Plan, the available capacity 
of these great laboratories is being turned over to work on the 
promising ideas of independent inventors everywhere, 


An Offer of Research Facilities 


to Inventive Americans Who Need Them 


The Sinclair Plan is opening up the Company’s great laboratories 
to every American who has an idea for a better petroleum product 


NVENTIVE Americans are often at a loss 
I today. Not because of any lack of ideas, 
but because of a need for expensive facilities 
to find out if and how their ideas work. 

This was no obstacle in our earlier days. 
The Wright Brothers designed their first air- 
plane with the help of a foot-square home- 
made “wind box’’—and the plane flew. 

In contrast, the man with a new idea in 
airplane design today often needs a super- 
sonic wind tunnel costing millions. 

In short, science and invention have be- 
come so complex that a man with an idea 
fora better product often needs the assistance 
of an army of specialists and millions worth 
of equipment to prove his idea has value. 

Within the petroleum field, the Sinclair 
Plan now offers to provide that assistance. 

Under this Plan, Sinclair is opening up 
its great research laboratories at Harvey, 
Illinois, to independent inventors who have 


sufficiently good ideas for better petroleum 
products or for new applications of petro- 
leum products. 

If you have an idea of this kind, you are 
invited to submit it to the Sinclair Research 
Laboratories, with the provision that each 
idea must first be protected, in your own 
interest, by a patent application, or a patent 


The inventor’s idea remains his own property 


If the directors of the laboratories select your 
idea for development, they will make, in 
most cases, a very simple arrangement with 
you: In return for the laboratories’ invest- 
ment of time, facilities, money and personnel, 
Sinclair will receive the privilege of using the 
idea for its own companies, free from royal- 
ties. This in no way hinders the inventor 
from selling his idea to any of the hundreds 
of other oil companies for whatever he can 
get. Under the Plan, Sinclair has no control 


over the inventor's sale of his idea to others, 
and has no participation in any of the 
inventor's profits through such dealings. 
Moreover, it is a competitive characteristic 
of the oil business that the new products 
adopted by one company are almost invar- 
iably adopted by the whole industry. This 
means that the very fact of his agreement 
with Sinclair should open up to the inventor 
commercial opportunities which might 
otherwise be hard to find. 

How to proceed: Instructions on how 
to submit ideas under the Sinclair Plan are 
contained in an Inventor's Booklet available 
on request. Write to: W. M. Flowers, Execu- 
tive Vice-President, Sinclair Research Labo- 
ratories, Inc., 630 Fifth Avenue, New York 
20, N. Y. for your copy. 

IMPORTANT: Please do not send in any ideas 
until you have sent for and received the 
instructions 


SINCLAIR — A Great Name in Oil 





—Contains strong, 
hair-like fibers which 
form an impervious 
mot through which 
mud will not poss 


—Especially prepared 
cellophone in floke 
form. 





—Contains thin flakes 
of carefully sized 
muscovite mica 


—A free-flowing 
blend of correctly 
sized Douglas Fir 
Fibers 


where you need these proved 


MAGCOBAR Products 


Magcobor hos developed these specific ‘formation plugging’ agents which 
effectively combat and correct conditions that contribute to loss of drilling fluids 
to the formation. In developing these products, Magcobar has given first considera 
tion to materials that will form an impervious bridge in fissures, cracks or holes, 
without fouling the valves and other working ports of the mud circulating pumps 
Magcobor also has available several plastic mud plugging agents for special 
problems. When you detect loss of circulation, coll your Magcobar field engineer 
for assistance and order out the Magcobar material of your choice from the 
Mogcobar dealer neorest you. There are now more than 300 strategically located 
dealers to serve you, night and doy with Magcobar’s line of complete Drilling 
Mud Service 


ALWAYS KEEP A FEW SACKS OF THESE MATERIALS AT THE RIG! IT'S GOOD 
INSURANCE AGAINST COSTLY DELAY AND HOLE HAZARDS 


MAGNET COVE BARIUM CORPORATION 
MALVERN, ARKANSAS . HOUSTON, TEXAS 


One of the Dresser Industries 
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Two Clark No. 2, 5-Impelier, Uncooled Centrifugal Compres- 
sors, eoch 750 bhp ot 8900 rpm, driven by a noncondensing 
steam turbine ot Sun Oi! Co.'s Marcus Hook, Pa. refinery, 
near Philadelphia 


UNIQUE APPLICATION BY SUN OIL’S 
MARCUS HOOK REFINERY 


V In feeding air system network 
Y In reducing steam pressure 


Wherever steam is available and space is at a 
premium, Clark Centrifugal Compressors have 
long been recognized as the ideal source of 
shop air. That’s a widely known fact! 

But, Sun Oi] Company’s refinery at Marcus 
Hook, Pa. had an additional and unique rea- 
son for recently installing two of these Clark 
units: They were most economically suited to 
the energy-conserving, cost-cutting 
which Sun had engineered, to accomplish in 


system, 


one continuous operation — (1) The feeding 
of the refinery air system network, and (2) 
The reduction of large volumes of steam from 


CLARK 


Centrifugal 
COMPRESSORS 


406 to 135 psi pressure. 

By using a noncondensing steam turbine to 
drive two Clark No. 2, 750 bhp (each), 5-Im- 
peller, Uncooled Centrifugal Compressors, 
and simultaneously reduce the steam pressure, 
Sun has come out far ahead. Available steam 
is being utilized to economically serve air re- 
quirements (the unit compresses 6150 cfm of 
air from the atmosphere to a discharge pres- 
sure of 110 psig). Large losses of energy, in- 
evitable with the use of a reducing valve, and 
high costs involved in the use of a boiler, have 
thus been avoided. 

To be current on the economy of Clark 
Centrifugal Compressors for your plant, write 
for literature and consult with your nearest 


Clark representative. 


SEE the difference in CLARK COMPRESSORS 
MIDGET ANGLE © RIGHT ANGLE © BIG ANGLE 
ELECTRIC-DRIVEN © CENTRIFUGAL 


CLARK BROS. CO., INC. 


One of the Dresser Industries 


OLEAN, N. Y. 


New York @ Tulse @ Houston @ Chicago @ Boston @ Washington @ Los Angeles @ Birmingham e Detroit « Salt Loke City @ Son Francisce 
london ¢ Poris @ Vorese, Italy @ Buenes Aires @ Coracas, Venezuela @ Lima, Peru @ Bogota, Colombic @ New Delhi, Indic 





W. realize that your choice of setting method is important, but before 
we consider the relative advantages of setting methods, let’s first 
analyze the end results obtainable with these tools after they are set. 

Among the important advantages of the Baker Retainer 
Production Packer are its uses... for production of one or more 
zones, simultaneously or separately, for acidizing, for re-pressuring— 
in fact for every application possible for a packer of its type. 

Features which appeal to veteran production men are its Com- 
plete Drillability—Resistance to High Temperatures—Corrosion- 
Resistant Construction— Freedom of the Tubing from the Packer— 
and Positive Anchoring against Upward or Downward Movement 


—Bottom-Hole Control with a flapper valve which closes when the 
BAKER tubing is removed, thereby eliminating necessity for killing the well 
with mud, under normal well conditions. 


MODEL ss The wide range of applications, permitting almost unlimited 


flexibility of subsequent operations, is based upon the construction 


RETAINER of the packer itself which employs two sets of opposing slips with a 


Buna N Rubber packing element between them. When it is set the 


PRODUC TION Packer virtually becomes a part of the casing, cannot be moved up 


or down, and prevents fluid or gas from passing between the packing 


PACKER element and the casing 


When the Setting Tool (used to set the Packer) is removed from 


PRODUC | the well, only the smooth-bore Packer remains, securely packed-off, 


= ready for selection and use of the proper Baker Packer Accessory 
NO 4}5-D ; Equipment to meet your production requirements. 

, f It is the extensive line of Accessory Equipment which permits 
installations to produce from one or more zones, as desired, through 
tubing or casing—and production changes often can be made with- 
out pulling the production string. Complete details of hook-ups, and 
illustrations of many unique, as well as routine installations, are 
contained in an 84-page Baker Packer Brochure, gladly sent to 
you upon request. 


SET ON TUBING OR A WIRE LINE 

In addition to conventional setting on tubing, the Baker Retainer 
Production Packer (Product No. 415-D) can be set on an electrical 
conductor cable by leading service organizations. See listing on 
Opposite page of the organizations in your area, ready to give you 
prompt, resultful service. 


BAKER OIL TOOLS, INC. 


Houston @ Los Angeles @ New York 
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When you select a Bridge Plug, it will pay you to choose the only one 
that is backed by more than 10,000 successful runs—at every depth 
—and for every known type of subsequent operation. The Baker 
Model “K” Retainer Bridge Plug ( Product No. 400) is an adaptation 
of the Baker Cement Retainer, which is too well known to all 
operators to require any explanation or comment. 

In addition to Anowing that you have selected the most successful 
Bridge Plug ever developed, you can choose the type exactly suited 
to your needs: The Cast Iron Type for permanent installations; the 
Magnesium Alloy Type for temporary service. 

Regardless of the size or type of Baker Bridge Plug selected, it 


ai 3 


will always drill out quickly and easily because all Baker Bridge Plugs 
are designed and constructed with “drillability” in mind. The very 
minimum cross-section of material need be drilled up, and even the 
slips are segmented to break up readily. 

You won't need to place cement above a Baker Bridge Plug to 
obtain a 100% leak-proof pack-off, but you'll be able to place 
cement there very simply if you desire with the Baker Model “B” 
Dump Bailer, which is run on a wire line at the same time as the 
Bridge Plug, when the wire line setting method is used. 

After the Baker Bridge Plug is set you can safely proceed with 
any future operations because it will withstand any pressures you 
dare impose upon the casing—remain permanently in position —and 
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drill out with no difficulty in record time —every time. 

CAN BE SET ON TUBING OR A WIRE LINE 

The Baker Retainer Bridge Plug can also be run and set on tubing- 
as thousands have been installed through the years—but to save 
wear and tear on tubular goods, as well as to save time on deeper 
settings, any of the leading service organizations listed below can run 
and set Baker Bridge Plugs on electrical conductor cables. Ask your 
service company for details—they will give you every cooperation. 


FOR IRE LINE SETTING CALL 


Byron Jackson Co. e Dowell, Inc. @ International Cementers, Inc. 
Lane-Wells Company e McCullough Tool Co. e Perforating Guns 
Atlas Corp. e Schlumberger Well Surveying Corp. e Welex Jet 
Services, Inc. e Well Perforators, Inc. e The Western Company 
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RETAINER 
BRIDGE PLUG 
PRODUCT 
NO. 400-D 
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DEEPER PENETRATION with Dowell Glass Gun in casing 


perforating or open hole shooting 


When you ask for Dowell Perfo-Jet Service, you receive 
the latest advancements in oil and gas well perforating! 
Now, Dowell has developed an expendable Glass Gun that 
allows the use of a larger explosive charge than do other 
types of guns. This extra detonating power gives you 
deeper, surer penetration of casing, cement and formation. 
Comparative “Standard drum tests” through 7’ O. D. casing 
into cement and sand show the following results: 


Hole Size 
5" diameter jet perforating gun 8%" ”" 
Dowell expendable Glass Gun 12” ar" 


Penetration 


DOwELL 


PERFO-JET SERVICE 


sk your nearest Dowell station for complete information on these Dowell services and products 
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and Bulk Inhibited Hydrochloric Acid. 


“First in Acidiz lg es 


. 8tnce 1932 





Dowell’s new Glass Gun is being used to good effect both 
in casing perforating and in open hole shooting. It consists 
of glass-enclosed, shaped charges in an expendable alumi- 
num carrier . .. and can be tailor-made to fit a wide 
variety of shooting patterns and hole diameters. This 
self-destroying gun assures you of a junk free hole. 


Call Dowell for your next perforating job. Ask for more 
detailed information on this new development. 


DOWELL INCORPORATED 
TULSA 3, OKLAHOMA 


Subsidiary of The Dow Chemical Company 





Loch 


dizing Service, Electric Pilot Services, Plastic Service, Chemical Scale Removal Service for 
at exchange equipment, Jelflake, Paraffin Solvents, Magnesium Anodes for Corrosion Control, 








FOR OIL INDUSTRY CHEMICAL SERVICE 





Offshore Production and 
Completion Practices 


~IN THE BAY MARCHAND 


AND MAIN PASS FIELDS 


By J. L. LeBlanc 
The California Co. 


Offshore Problems Do Not End With the Exploratory Drilling Phase. Production and Transportation of 
Crude Oil From the Well Site Offers Additional Mountainous Problems. Changes in Completion Practices 
Must Be Effected in Some Cases. These and Additional Problems Are Discussed Here for Two Fields Off 


the New Orleans Coastline. 


ABSTRACT 


While the tremendous job of initiating exploratory drilling 
in the Gulf of Mexico was in progress, thought was also being 
given to the problems which might be encountered after pro- 
duction was discovered. Those working on these problems soon 
began to realize that the production and transportation of 
crude oil in the offshore areas offered as many difficulties as 
the drilling. As The California Co. discovered fields on the 
second and fourth offshore wiidcats, production problems soon 
became reality. Concentration of drilling in these two fields, 
the Bay Marchand and the Main Pass Fields, accentuated 
production problems which seemed to have no possible eco 
nomic solution. Fig. 1 indicates the general location of the 
two fields. The Bay Marchand Field is located approximately 
65 miles south of New Orleans or about 3 miles offshore from 
the mouth of Bayou Lafourche, whereas the Main Pass Field 
is located approximately 75 miles southeast of New Orleans 
on the east side of the Mississippi River Delta and approxi 
mately 3 miles offshore. 

In conjunction with the surface problems, subsurface prob 
lems began to arise from wells directionally drilled to shallow 
unconsolidated sands. Completion practices had to be changed 
to meet the new situation. A discussion of these subsurface 
problems and practices will be given in the latter part of 
the paper 


SURFACE PRODUCTION EQUIPMENT 
AND PRACTICES 


Many schemes had been devised for the separation and 
transportation of offshore oil. All of these schemes seemed to 
be physically possible but few were even close to being eco 
nomically practical. 

Although these plans had been given careful consideration 
nothing concrete had been agreed upon even up to the time 
the discovery well in the Main Pass Field was scheduled to 
be brought in. A temporary means of testing this well con- 
sisted of a separator installed on the drilling platform, a 
65-bbl standby fresh water tank as a gauging tank and a 
barge moored alongside. These temporary installations served 
until drilling on the fourth well had begun. With the barge 
moored alongside, the drilling structure was continually being 
damaged. An anchoring system was devised for the barge, 
but this also proved impractical. Production subsequently dis 
covered in the Bay Marchand Field was being handled in the 
same manner except with even greater difficulties. One can 
well imagine the hazards and difficulties encountered with such 


Paper presented at the Petroleum Branch Fall Meeting in New Orlean 
La., Oct. 4-6, 1950 
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a production system in the Bay Marchand Field where 4 ft 
to 6 ft waves are common and 10 ft and 12 ft and even higher 
waves are not uncommon. We soon began to see that if pro- 
duction was to be placed on a daily basis a new system had to 
be inaugurated. Shallow producing horizons with their attend- 
ant sand problems further emphasized the necessity of aban- 
doning the productive systems in use. Sufficient wells had been 
drilled or were in sight to justify economically the installation 
of 4in. pipe lines to shore in both the Main Pass and Bay 
Marchand Fields. Pipe lines to shore based tankage and treat- 
ing facilities had been judged the most economical and reli- 
able. Accordingly, two 250-bbl tanks, separators and pumps 
were installed on the Main Pass “A” platform and a 4-in. pipe 
line was installed from this structure to a shore base where 
two 10,000-bbl tanks were erected. In the Bay Marchand Field 
flow lines were installed from the “B” and “C” structures to 
the “A” structure where two 250-bbl tanks, separators and 
pumps were erected and a 4-in. line connected these facilities 
to a shore base where four 10,000-bb] tanks were located. 
These facilities soon proved to be well worth their installa- 
tion costs. Hazards to field personnel were cut tremendously 
Production was placed on a daily basis and allowable produc 
tion was being made for the first time. These facilities con 
and were thought to be sufh 
until the hurricane 1949. It 
that they embodied many drawbacks in 
addition to some uneconomical aspects already foreseen. When 


tinued to serve adequately 


ciently practical season of 


was then realized 


FIG. 1 — MAP SHOWING LOCATION OF BAY MARCHAND AND MAIN 
PASS FIELDS 
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FIG. 2 — WATER KNOCKOUT SEPARATOR. 

a hurricane originates in the Gulf of Mexico, waves build up 
to an alarming height in the Bay Marchand Field, even before 
the hurricane has been located. By the time hurricane notices 
are received the waves are of such height (15 fi and higher) 
that it is an impossibility to get personnel on and off the 
various production platforms. Wells must be left open with 
reliance being placed on flow line safety valves and tubing 
storm chokes, neither of which is effective. In addition to these 
features of insecurity, the system offered the possibility of 
becoming too cumbersome and requiring excessive operating 
manpower as development progressed. 

It was agreed that a pipe line to shore was the best means 
of transportation but inasmuch as low tubing pressures and 
sanding conditions precluded the use of well energy it was 
necessary to have separating and pumping equipment offshore. 
Therefore work was started on the design of equipment which 
would meet all of the usual production requirements, in addi- 
tion to being automatic to the extent that manpower could be 
reduced extensively. This work resulted in providing the sys- 
fems now in use. 

Fig. 2 provides a cross section view of the horizontal sepa- 
rater which is the piece of equipment around which all of 
the system is built. This separator is basically a pressure vessel 
5 ft in diameter and 20 ft long, having a volume of approxi- 
mately 70 bbl and constructed to withstand 125 psi operating 
pressure. In addition to being a separator, it acts as a surge 
tank, test tank and free water knockout. A complete produc- 
tion unit for a platform consists of two of these horizontal low 
pressure separators (one for bulk and one for test), two con- 


ventional horizontal high pressure separators, a pipe line 
pump. well head safety valves and necessary pneumatic 


controls, 

Fig. 4 is a view of a complete unit in operation on the Bay 
Marchand “C™ 
the platform is seen in the foreground. Each low pressure 
a high level 


platform. One of the five producing wells on 


separator is equipped with two safety features: 
controller which shuts in the wells in the event that the oil 
level in the separator rises above normal; and a high pressure 
controller which shuts in the wells in the event the pressure 
within the separator becomes too high. Each high pressure 
separator is equipped with a high level controller which shuts 
in the wells in the event the oil level in the separator becomes 
too high. High pressure controllers have not been installed on 
the high pressure separators because of the separator capacity 
to withstand well head pressures. Pilots on these controllers 
pneumatically -activate piston operated Christmas tree safety 
valves which shut in the wells in the case of pump or other 
failures. The float on the low 
controls which start and stop the pipe line pump at the appro- 
priate time. Electrically driven pumps are used in the Bay 
Marchand Field where electrical 
through submarine cables. In the Main Pass Field, pumps are 


pressure separator operates 


power was made available 


driven by pneumatically controlled internal combustion en- 


gines. 
> 


Periodic well tests can now be made with ease as the LI 


test separator has sufficient volume to permit 2 or 3-hour 
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testing of wells. It is no longer necessary to provide gaugers 


and living quarters on each structure as was the case under 
the old systems. Two gaugers are now employed in each field 
whereas six would have been required under the old system. 
The cost of maintaining standby boats has also been reduced 
proportionally which makes a substantial reduction in daily 
operating costs. These automatic producing facilities have so 
far been practically trouble free, and their operation has been 


highly pleasing. 


SUBSURFACE PRODUCTION AND 
COMPLETION PRACTICES 

Bay Marchand Field 

The Bay Marchand Field is situated on the largest known 
piercement type salt dome in the world. The salt has pene- 
trated to very shallow depths, being contacted at 2,110 ft in 
the first All production to date been found in sands 
overlying the dome. No deep flank wells have been drilled as 
yet except on the northern flank which extends inshore and 


well. 


can be drilled by standard inland drilling barges. The northern 
flank has not been productive. Numerous tangential and radial 
faults complicate the subsurface geology. This in conjunction 
with slant holes makes the area a geologist’s nightmare. 

To date, six drilling platforms have been constructed in 
the field in an average water depth of 30 ft. From these plat- 
forms 27 wells have been completed as producers giving a 
total daily production of approximately 6,200 bbl. Production 
is from 13 different sands ranging in depth from 2.350 ft to 
5.700 ft. The reservoirs being produced are of the water drive 
type with reservoir pressures being normal or near normal. 
The oil gravity ranges from 21° API to 41° API. Sixteen of 
the wells are producing by natural flow and four are produc- 
ing by gas lift. Fourteen of the wells have been gravel packed 
to eliminate sanding troubles. 

Two of the six platforms in the field were constructed so 
that two drilling rigs could be placed on the platform and 
both rigs could operate simultaneously. As all of the wells 
except one on each platform were drilled directionally. it was 
of major importance to plan development drilling carefully in 
order to prevent wells from running into each other. Fig. 3 
indicates that even with careful planning two such situations 
came very close to arising on the “D” platform. Each circle 
on the sketch represents a vertical depth of approximately 
100 ft. As the wells were drilled on 8 ft centers. 


was 


every effort 


made to keep the vertical sections of the hole straight. 
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FIG. 3 — PLAN VIEW, ‘’D’’ WELLS, BAY MARCHAND FIELD. 
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but as Fig. 3 shows the holes wandered considerably in spite 
of these efforts. Fig. 3 shows surface projections of the upper 
ections of the holes only as this is the only section where 
danger of collision existed. 


Development in the Bay Marchand Field has been on an 
approx.mate 20-acre spacing program. Horizontal displace 
ments of the wells have therefore averaged between 900 ft and 
1,800 ft. If possible the wells are drilled at angles somewhere 
between 18° and 25° off vertical as this range is considered 
to be the best for directional control, but on occasions wells 
were drilled at angles of 30° and 35° off vertical and one one 
wildcat in the field an angle of 54° was reached without undue 
difficulty. The use of oil emulsion mud to reduce torque and 
drag in these directional holes has now become a standard 
practice. From 10 per cent to 15 per cent by volume of diesel 
oil is added to the mud. No particular effect has been noticed 
on sidewall samples and gas and water contacts are as readily 
established as they are in holes drilled without oil emulsion 
muds. 

The unconsolidated nature of the shallow sands encountered 
in most cases demanded that some means of sand control be 
established during the completion programs. The use of gravel 
pack was considered the best means of sand control available 
for the conditions present in the Bay Marchand Field. Three 
types of gravel pack jobs have been employed in the field 
(1) the single stage or inside casing pack, (2) the two stage 
pack inside casing where on the first stage gravel is squeezed 
out through the perforations and on the second stage the well 
is packed inside the casing, and (3) the open hole pack where 
a section of the casing is milled out, the hole underreamed 
and an open hole pack performed. All three types of jobs 
have proved successful from the standpoint of sand control but 
have varied widely from the standpoint of well productivity 
Productivity indices have been established on all of the packed 
wells and vary from an average of .350 for the single stage 
packs to an average of 1.40 for the two stage packs and 4.3 
for an open hole pack. This compares with an average PI of 
6.0 for the unpacked wells in the field. From a productivity 
standpoint the open hole pack is definitely the best, but the 
costs of milling out casing and underreaming has until recently 
discouraged this approach when the other methods would suf 
fice. Completions involving serious sand problems have until 
recently been standardized on the two stage pack, which is a 
compromise between cost and best results. However, three 
more open-hole packs, in addition to the above mentioned 
gravel packs. involving mill-out of casing and underreaming 
of sand have just been made with the anticipated good pro 
ductivity and solution of sand problems and at a cost that, 
although high, is justifiable. The possibility of setting casing 
on top of a sand or the use of aluminum alloy casing for 
easier milling thus making open hole packs cheaper have both 
been given consideration, but the particular well features 
necessary for use of either procedures have not as yet been 
present. 

One of the detrimental characteristics of a single stage pack 
is that it will gradually plug up. reducing flowing tubing pres 
sures and production in proportion. This has been particu 
larly true in the Bay Marchand Field where flowing tubing 
pressures, even though still reasonably high, are declining at 
a faster rate than anticipated. The one open hole pack per 
formed in the field has as yet given no indication of declining 
flowing tubing pressures or production even though approxi 
mately 150.000 bbl have been produced from the well. This 
compares with the single stage packs where initial declines 
been noted within the first 20,000 bbl of cumu 
lative production, Most of the successful two stage jobs per 
formed in this field and in the Main Pass Field have been 


have usually 
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FIG. 4— COMPLETE UNIT IN OPERATION, BAY MARCHAND “C” PLAT- 
FORM 


fairly recent, and no concrete information has been gathered 
as yet, although it is expected that these packs will have more 
favorable production histories than the single stage packs. 
\ change in the size of gravel used has given no noticeable 
change in production characteristics. Gravel ranging in size 
from .060 to .090 in. was used on the first inside packs per- 
formed. This size is still being used on the first stage of the 
two stage packs but .040 x .060 in. gravel is now used inside 
the casing. The possibility of using even smaller gravel has 
been given consideration and will probably be tried in the 
future. 

On several occasions diesel oil has been pumped down the 
tubing and through the gravel pack in an attempt to back 
wash the pack and thus eliminate the clogging effect. In gen- 
eral, if pumping pressures are not excessive, about 50 bbl of 
diesel are pumped through the pack and back into the forma- 
tion. Results have not been phenomenal but have been satis- 
factory to the extent that the expenses involved are considered 
well worthwhile. On two occasions flowing tubing pressures 
have been restored to normal. The results no doubt would be 
if selective washing could be done 
but this would require the use of a rig. 

Main Pass Field 

Subsurface production practices in the Main Pass Field 
have been substantially the same as those employed in the 
Bay Marchand Field and for that reason will not be discussed 


much more satisfactory 


in detail. Three drilling platforms have been constructed in 
the field from which 1] wells have been drilled? Nine of these 
wells have been completed as producers in four sands ranging 
in depth from 5,500 to 8.150 ft. Daily production is approxi- 
mately 3,000 bbl. 

The field is situated on a deep-seated salt dome structure. 
All reservoirs being produced are of the water drive type with 
reservoir pressures and temperatures being normal. Faulting 
is general in the area but not to the extent found in the Bay 
Marchand Field. 

Sand problems are present but have not in general been as 
prevalent as in the Bay Marchand Field. Three of the eight 
producers have been completed with gravel pack. Of these 
three packs two are of the single stage inside casing type and 
the third is the two stage type. Productivity indices of the 
packed wells are very similar to those found at Bay Marchand 
Only a slight decline in tubing pressures of the packed well. 
has been noticed but these pac ks have been too recent to war- 
rant any conclusions. Packing procedures have been identical 
to those Marchand so no wide differences in 

es & 2 


used at Bay 


results are expected, 
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ENGINEERING DEVELOPMENTS IN LOCAL SECTIONS 


Members of the Student Chapter at 
the University of Southern California 
met last month for their annual all-day 
meeting and heard, in addition to a 
brace of technical papers by graduates 
and undergraduates, a sharp criticism 
of the present Administration’s eco- 
nomic policies. They were told that the 
achievement of any great financial suc- 
cess under present economic policies is 
“more or less impossible.” Author of 
these criticisms was W. H. Geis, con- 
sultant, and banquet speaker for the 
occasion, 

Earlier in the day, awards were made 
by the Pacific Petroleum Chapter and 
the AIME Petroleum Branch office in 
Dallas to first and second place win- 
ners in the graduate and undergraduate 
paper contest. In the graduate petro- 
leum engineering division, Barry L. 
Evans took first place and George 
Ichillingarian took second place. In 
the undergraduate division, R. C. Kirch- 
man was awarded first place. The Pa- 
cific Petroleum Chapter gave all win- 
ners a slide rule and the Branch re 


warded first place winners with copies 





A USC Student Chapter Report 


of Subsurface Geologic Methods. Chair- 
man of the Awards Committee was Paul 
Andrews of Signa! Oil and Gas Co 

At the evening banquet, Geis glumly 
warned the young engineers that they 
“probably will suffer be greatest frus- 
trations if the present controlled econ- 
omy, which in fact is a form of eco- 
nomic socialism, continues.” He com 
pared the economic programs of Social- 
istic England and the United States in 
this manner: 

So far England has nationalized, or 
rather expropriated from private enter- 
prise, the coal, steel, railroad, banking, 
and some other industries, including the 
general direction of labor as a form of 
industry. That means that the English 
government receives all the net revenue 
or the net profits from 50 per cent of 
its industries. This year in the United 
States the government, through the oper- 
ation of the federal income tax laws, 
will receive not less than 55 per cent, 
and in some cases as much as 77 per 
cent, of the net revenue, that is, the net 
profit, from 100 per cent of the United 
States’ industries. In comparative terms, 


USC STUDENT PROGRAM — Siudent paper contests held last month at the 
University of Southern California produced three winners in the petroleum engi- 
neering division. Shown above are, L to R, Barry L. Evans, first place winner in 
the graduate class; R. C. Kirchman, first place winner in the undergraduate 
class; Paul Andrews of Signal Oil and Gas Co., Awards Committee chairman: 
W. H. Geis, banquet speaker; and George Tchillingarian, second place winner 


in the graduate class. 


1% . . . SECTION 1 


JOURNAL OF PETROLEUM TECHNOLOGY 





Some Economic Aspects of the Oil Industry 


while 50 per cent of English industry 
works 100 per cent of the time for Eng- 
land, 100 per cent of Americaa industry 
works 55 per cent, or 77 per cent, of the 
time for the American government. Un- 
der such a comparison, which govern- 
ment has turned further towards so- 
cialism? 

In addition, another avowed purpose 
of the socialists and communists is to 
destroy private capital. Mexico, at least, 
made no bones about it. Mexico simply 
expropriated all the properties and as- 
sets of the oil companies and so did 
Red China with British Asiatic Oil Co. 
These were direct and positive actions 
But the more favored way to destroy 
capital is by the slower, but just as 
sure, method of levying a tax on vapi 
tal. It is a well known economic prin 
ciple that whenever capital, as such, 
is taxed it must eventually follow that 
the capital is destroyed. 

In this regard Britain is smart, for 
with all of its advertised socialistic 
trends it has never yet levied a tax on 
capital or on capital gains; nor has it 
ever levied a tax on dividends, which 
are a reward for the use of capital by 
industry. 

On the other hand, the United States. 
by indirection and by various devious 
means, levies a tax on 85 per cent of the 
capital in the mining industry and on 
72% per cent of the capital of the oil 
industry. The government does this con 
trary to basic, sound economic prin- 
ciples, simply by classifying the return 
of capital as ordinary income. As a 
result, the independent mining indus- 
try which flourished so well in its early 
history is practically non-existent in 
the United States today because its 
capital is gone and investors will not 
risk capital against such odds. It is, of 
course, true that a very few powerful 
mining companies remain which, largely 
because of their size, have been able to 
revive their capital through outside 
investments. In addition, an income tax 
is levied against all dividends in the 
United States. and the Treasury De 
partment has now asked Congress to 


Continued on Page 2, Section 2 
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History and Development 


of the AAODC Training Programs 


The most pressing need in the drill- 
ing industry today is for trained per- 
sonnel, believes A. W. Thompson of 
Thompson-Carr, Inc., Houston, presi- 
dent of the American Association of 
Oilwell Drilling Contractors. Thompson 
outlined the development of the AAODC 
training programs before a meeting of 
the Permian Basin Section May 21. 

Training is reflected, Thompson said, 
in ways other than the technological 
end of the work. There is also developed 
an improved morale and attitude. 
greater teamwork among the men and 
an intangible quality of workmanship 
that leads to more and better footage of 
hole per dollar spent. Thompson attrib- 
utes the increase in drilling efficiency in 
recent years partly to the association’s 
research and educational programs, all 
of which were begun in the Permian 
Basin area. 

Drill pipe corrosion was the first 
problem to be approached by the con- 
tractors. Battelle Memorial Institute, 
Columbus, Ohio, was selected to con- 
duct a program of research, which was 
financed by individual subscription of 
contractors. A Research Committee of 
AAODC was appointed to work with 
Battelle, and the program was initiated 
in July, 1945. 

One of the conclusions of the studies 
was that corrosion of drill pipe was 
closely related to the muds used. In 
September, 1945, John Woodruff, State 
Coordinator of Petroleum Industry 
Training at the University of Texas, 
proposed to the association a mud re- 
search and mud training program in 
collaboration with the University. This 
program went into effect in February, 
1946, and in April of that year the as- 
sociation made a grant of $6,000 to the 
University of Texas for expansion of 
the program. 

Later in the year the first Research 
and Education Committee was formed 
With the increase of deep drilling, in- 
dustry began to feel the need for trained 
men more acutely. In order to give ba- 


sic training to young men entering the 
business as drillers and toolpushers, a 
Drilling Practices Course was organ 
ized in February, 1947. The first course 
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Reported by Brad Mills 
Permian Basin Section 


was set up at Kilgore College in co- 
operation with the University of Texas. 
A complete drilling rig was donated for 
this work by various drilling contractors 
and supply firms. During that year the 
present programs were expanded to in- 
clude further grants to the University 
of Texas and additional Drilling Mud 
Control Schools in various parts of the 
country. 


Two new training programs were or 
ganized in 1948. A supervisory pro- 
gram, “Accident Prevention Through 
Better Supervision,” was established on 
a conference basis. a nd conferen¢ es 
were held in Texas. Wyoming and Ok- 
lahoma during that year. A series of 
Management Institutes was organized 
to begin in 1949, covering functions of 
management, personnel problems, acci- 
dent prevention, cost factors, security 





RUPERT C. CRAZE, research special- 
ist for Humble Oil and Refining Co., 
spoke on “The Performance of Lime- 
stone Reservoirs” before the monthly 
meeting of the Gulf Coast Section 
June 18. He discussed the fundamental 
principles governing oil and gas pro- 
duction from sandstone and limestone, 
and the applicztion of analytical tech- 
niques used for sandstone reservoir 
studies to the study of limestone res- 
ervoir performance. 
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for the drilling organization and legal 
problems of the drilling contractor. A 
Drilling course was also 
planned. 


Engines 


The first Management Institutes in 
1949 had an advance registration of 
nearly 500 men. The Mud Schools and 
the Supervisory Training Programs were 
continued at the same time. Assistance 
from several insurance companies was 
secured for the Supervisory Training 
Programs. Supply and manufacturing 
companies donated a complete engines 
laboratory for the first Drilling En- 
gines Course, organized in September, 
1949. A new Educational Administra- 
tive Committee had been appointed by 
this time to coordinate the expanded 
programs and supervise expenditures. 


During 1950, efforts were made to 
expand the Employer Unit Plan of 
Supervisory Training, which had been 
started in the Permian Basin with the 
Carl B. King and W. H. Black Drilling 
Companies. Management Institutes were 
held in 11 different drilling centers for 
215 persons. Topics for this second se- 
ries, selected from a poll taken at the 
first institutes, were contractor-company 
relationships and insurance problems. 

During 195] the Mud Schools have 
been expanded, the Supervisory Train- 
ing Program has been enlarged to in- 
clude California and the Rocky Moun- 
tains, and Oklahoma units have in- 
creased. The Gulf Coast unit has been 
added. A new course to supplement the 
Employer Unit Plan is known as Field 
Institutes. It is a specialized training 
aimed at toolpushers and drilling su- 
perintendents, providing conferences on 
pecific drilling problems. A new proj- 
ect is underway to publish an illus- 
trated booklet on oilwell drilling termi- 
nology, written primarily for office per- 
sonnel in drilling and oil companies. 

Currently, Mud Schools have been held 
in Kansas, Casper, Wyo., and West 
Texas; Field Institutes are being held 
once a month in the Permian Basin and 
in Casper; and Employer Units are 
operating in the Permian Basin, Gulf 
Coast, Oklahoma, California and the 


Rocky Mountain areas. R. F2 
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Predicting Reservoir Performance — II 


Reported by John Clayton Jacobs, Jr. 


The two final Spring meetings of the 
study group of the North Texas Section 
were held on May 9 and May 23, led 
by A. B. Dyes of the Research and 
Development Department of the Atlan- 
tic Refining Co. These two discussions 
dealt with the application of material 
balances to the prediction of reservoir 
performance. Primarily emphasized was 
the means of evaluating the magnitude 
of a reservoir water drive and the man- 
ner in which flash and differential P-V 
data are employed in material bal- 
ances. 

In developing the optimum method 
of producing an oil or gas reservoir, the 
magnitude of any water drive which 
may be present and its effectiveness in 
the future must be known. The water 
drive at any given date can be evaluated 
by material balance if the quantity of 
in place material, production and pres- 
sure histories, and fluid properties are 
accurately known. However, since the 
volumetric estimate of in place material 
is only approximate, the definition of 
the water drive by this means. if done 
early in the life of a field, will be in- 
accurate and can lead to large errors 
in the prediction of the future effective- 
ness of the water drive. The uve of 
many material balances employing sev- 
eral different degrees of water drive 
with several values of in place material 
to define the pressure decline will per- 
mit the selection of the in place mate- 
rial and water drive to yield a minimum 
deviation of the calculated pressure de- 
cline from that observed in the field. 
These quantities can then be used in 
caleulating the future behavior of the 
reservoir. 

The future water drive cannot be 
obtained accurately by an extrapola- 
tion of the coefficient in the Schilthuis 
or “short form” Hurst equation except 
under very special reservoir conditions. 
These conditions are not generally pres- 
ent until late in the producing life of 
the reservoir. In predicting the amount 
of water influx with time. either the 
electric analyzer or the “long form” 
calculation procedure developed — by 
Hurst should be employed to give the 
most reliable results. 

In predicting the performance of a 
reservoir, utilizing material balances. 
the P-V properties of the reservoir fluids 
must be known. These are described by 
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North Texas Study Group 


two types of P-V data, one obtained by 
differential liberation of gas from the 
reservoir fluid over small pressure 
changes at reserveir temperature and 
the other by flash liberation of gas from 
the reservoir fluid by one reduction in 
pressure to the separator pressure and 
temperature as present in the field. Both 
differential and flash data must be em- 
ployed in making balance 
studies. The stock tank oil and separa- 


material 


tor gas production as reported from the 
field should be converted from flash to 


differential production and employed in 
a material balance using differential 
P-V data to describe the performance 
in the reservoir. Likewise, predictions 
in the future would be obtained from 
material balances employing differential 
data by converting the predicted differ- 
ential production from the reservoir to 
flash production in the stock tank. The 
use of only one type of P-V data can 
lead to appreciable errors in the esti- 
mation of in place material and reser- 


voir performance. x * * 


Making Research Results Useful 


Reported by Herman Dykstra 


Pacific Technology Group 


The fourth lecture of the present se 
ries on research being conducted by the 
Pacific Technology Group was given by 
Howard C. 
on May 28, at the General Petroleum 
Auditorium in Los Angeles. 

Pyle began by defining research ac- 


Vie, petroleum engineer, 
Pyle, petrol 


cording to Webster and by classifying 
research as pure basic, or long term, 
research and industrial, or short term, 
research. Long term research is mainly 
concerned with the investigation of 
fundamentals, whereas short term re- 
search is mainly concerned with an 
swering specific problems as are often 


found in industry. An example of a 


long term research is API Project 37, 
“Fundamentals of Hydrocarbon Beha 
vior.” Industry does not usually take 
such a research project because the 
application of results is not immediately 
apparent. 

As an example of short term research 
Pyle described the development of mud 
scratchers. Someone suggested the use 
of a mud scratcher and accordingly re- 
search was carried out to determine the 
effectiveness of a mud secratcher. Within 
a few weeks results were obtained which 
showed that mud scratchers were feas- 
ible. 

Actually many research projects are 
both long term and short term projects 
Useful results are obtained during the 
research work but, nevertheless, the 
project is extended for years. 

All research eventually pays off and 
has a useful objective. Research is fi 


nanced by industry, and employers want 
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results and facts, and many times may 
want to expand on the basis of results 
already obtained. 

Pyle concluded by giving several rules 
that he believed researchers should fol- 
low in order to get maximum usefulness 
from research results. 

1. Researchers should be enthusiastic 
about their work and should keep their 
research objectives foremost in their 
conversation. They should avoid the 
temptation of talking about interesting 
phases of their research that are not 
directly related to the objectives. 

2. Researchers should speak in sim- 
ple terms easily understood by man- 
agement and those who hold the purse 
strings. 

3. Researchers should appreciate eco- 
nomics of research and should know the 
cost of carrying on research. On the 
should not 
measure all research in terms of dol- 


other hand management 
lars and cents. 

}. Research reports should be writ- 
ten for the non-researcher in simple 
straightforward Technical 
supporting results should be placed in 
an appendix. This includes all differ- 


language. 


ential and integral signs. 

5. Use of graphs and illustrations are 
helpful to bring out facts and figures. 

6. The researcher should keep man- 
agement informed of results. 

7. Research men should be prompt 
in presenting their ideas and results of 
their work. Research is big business 
and should be handled as such. * * * 
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THE QUANTITATIVE APPLICATION 
OF RADIOACTIVITY LOGS 


ROBERT E. BUSH, LANE-WELLS CO., HOUSTON, TEX., AND E. S. MARDOCK, MEMBER AIME, WELL SURVEYS, INC., 
TULSA, OKLA. 


ABSTRACT 


The object of this paper is to report recent developments 
in the quantitative interpretation of radioactivity logs. The 
use of reference lines is described in the application of the 
new zero radioactivity reference line for the determination of 
neutron derived porosities is discussed. Illustrated examples 
are shown of logs from the reef limestone of the North Snyder 
Field and from the Lansing-Kansas City limestone of central 
Kansas. 

\ method of using an interrelation factor calculated from 
the neutron and gamma ray curves, called the neutron pro 
ductivity number, for the estimation of productivity 
sandstone formations is discussed. When used in a relative 
manner the neutron productivity number has aided in pre 


from 


dicting whether a given zone will produce gas or liquid or 
whether it is too shaly to produce. Examples of logs from 
wells in the Freites District. Venezuela. and in East Texa 
are shown to illustrate the application of the neutron produ: 
tivity number. 


INTRODUCTION 


The petroleum industry was introduced to some of the pre 
liminary results of investigations into the quantitative inter 
pretation of radioactivity logs in October, 1949." It was stated 
at that time that the methods presented were entirely relative 
because they involved the use of reference points obtained 
from the radioactivity log, i.e., shale reference lines, dens 
zone reference lines, or a combination of both shale and dense 
zone reference lines. These relative methods offered many ad 
vantages to the industry in interpreting radioactivity logs in 
terms of porosity with particular application to limestone and 
dolomite reservoirs. However, these relative methods did not 
offer the wide range of applicability desired because the e-tab 
lishment of reliable reference points from shale in some areas 
were difficult. It was recognized that a more reliable means of 
establishing reference lines was vitally needed. Such a method 


‘References given at end of paper. 

Manuscript received in the office of the Petroleum Branch Sept. 8 
Paper presented at the Fall Meeting of the Petroleum Branch 
Orleans, La., Oct. 4-6, 1950 
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is now available through the introduction of an instrumental 


zero device. It is the purpose of this paper to discuss the 


application of this new zero reference to the quantitative 
nterpretation of radioactivity logs. 

Previously all applications of quantitative interpretations 
of radioactivity logs have been made in limestone reservoirs 
with little or no emphasis placed on sandstone reservoirs. A 
brief resumé of investigations into the quantitative interpreta- 
tion of radioactivity logs in sandstone reservoirs is presented 
here. This discussion deals primarily with the interrelation of 
the gamma ray and neutron curves and the response of the 
latter to the density of the fluids filling the pore space in 
relatively homogeneous sandstones of uniform porosity. The 
authors believe that the methods presented herein are only 
applicable locally to specific fields and reservoirs, and that 
they probably do not offer a wide range of applicability. 

However, the authors believe that the wider use of the new 
instrumental zero references in sand areas will reveal more 
information about the specific response of the gamma ray and 
neutron curves to changes in lithology and saturation. 

Two applications of quantitative interpretation using the 
new instrumental zero references are presented. Examples are 
given which show the comparison of logs and porosity in the 
limestones and dolomites of West Texas and Kansas. Logs of 
wells in Venezuela and East Texas are compared with the drill 
stem test data, production data, ete.. to illustrate the inter- 
relation of the gamma tay and neutron curves and the use of 
interrelations to determine certain 


methods involving these 


saturation variables. 


THE QUANTITATIVE APPLICATION OF 
RADIOACTIVITY LOGS TO LIMESTONE AND 
DOLOMITE RESERVOIRS 


rhe theoretical basis for the gamma ray and neutron curves 
has been discussed elsewhere.” The method of recording the 
ogs and apparatus involved has also been thoroughly covered 


ind will not be 


disc ussed bere. 


It has been shown that the neutron curve responds expo- 


nentially to liquid filled porosity over the 11% to 35 per cent 
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GAMMA RAY 
SENSITIVITY 6.32" 


NEUTRON 
SENSITIVITY 4.13” 


INSTRUMENT ZERO LINE 
SHALE LINE 
SHALE LINE 











FIG. 1 — RELATIONSHIP BETWEEN RELATIVE REFERENCE LINES AND 
THE ZERO RADIOACTIVITY REFERENCE LINES FOR 7” CASED HOLE 


range in limestones and dolomites.’ It is assumed therefore 
that the relationship between the neutron curve deflections and 
porosity for a given zone, can be determined graphically by 
using this relationship to weight the average line drawn 
through a plot of the scattered values of the neutron curve 


deflections and the corresponding porosities. 


Introduction of Instrumental Zero 
Reference Lines 


Previous work has made use of relative methods in which 
from the shale line of the 
neutron curve. However, shale lines were not always easily 
established and _ it 
method of obtaining a better zero reference was needed. A 


reference lines were established 


was therefore quite apparent that some 


zero locating device has now been added to the existing logging 
equipment. It is now possible to determine the position of zero 
radioactivity on the log at any time during the logging opera 
tion. Fig. 1 shows the relationship between the instrument 
zero and the shale reference lines for the gamma ray and 
neutron curves in a_ typical limestone area. To acquire more 
adequate control of the sensitivity, a sensitivity selector is 
now available which permits accurate measurements of the 
sensitivity and the setting of any desired sensitivity. The addi 
tion of this equipment makes it possible to standardize sensi 
tivities for given logging conditions: i.e., under given casing 
and core hole sizes. 

To illustrate the application of quantitative methods to lime 
North 
Tex.. was chosen because core 
same well for both cased hole 
2 shows the correlation of the 


stone reservoirs, a well in the reef limestone of the 


Snyder Field, Scurry County. 
data were available within the 
and open hole conditions. Fig. 
radioactivity log and the core porosity log of a standard Oil 
Co. of Texas well in the North Snyder Field. 
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Preparation of the Radioactivity Log 
For Quantitative Interpretation 


le prepare the 
tion the 


radioactivity log for quantitative interpreta 
location of the absolute zero reference line is deter- 
mined at one or more points within the zone of interest and 
throughout the remainder of the log on beth the gamma ray 
and neutron curves. The porosity determined by core analysis 
is then plotted on the same vertical scale as the radioactivity 
log and correlated core by core with the neutron curve. Where 
physical bore hole conditions such as casing or bit size are 
changed within a well it is 


necessary to prepare a separate 


quantitative interpretation chart for each section of the well. 
\ data sheet showing the depth. the gamma ray curve deflec- 
tion from zero, the neutron curve deflection from zero, and 
the porosity of the corresponding core sample is then prepared, 
throughout the cored section, until all of the data necessary to 


prepare the neutron derived porosity curve are recorded 


Neutron Derived Porosity Curve 


assuming that the neutron curve responds in an expo- 
nential manner to liquid filled porosity. the establishment of 
facilitated if the data 
ere plotted on semi-log graph paper with the neutron curve 


the neutron derived porosity curve is 


deflection as the abscissa (linear scale) and the per cent core 


porosity as the ordinate (log seale). The neutron derived 


porosity curve is determined by line 


that 


drawing the straight 
When the 
distance of pen 


best fits the scattered points. neutron curve 


travel 
essary to adjust for changes of sensitivity 


deflections are plotted in units of 
meches tf is nec 
when applying the data of one log to that of another log. The 
to the deflection and 


therefore the neutron derived porosity curve may be adjusted 


relationship of the sensitivity is linear 


to another sensitivity by multiplying it by the ratio of the 


new to the old sensitivity. 


However, if the neutron curve deflection is divided by its 


sensitivity, an invariant quantity is obtained which is applic- 
This quantity has been called a standard unit, 
be plotted in 


either standard units or inches of deflection as desired. Figs. 2 


ible to any log 


and the neutron derived porosity curve may 


and 3 show an example of the use of standard units, while an 
f xampl of a neutron derived porosity curve plotted in inches 
Figs. 5 and 6. The value in standard 


units for any neutron curve deflection may be found by divid- 


of deflection is shown by 


the curve deflection (in inches) at the given point by the 
sensitivity of the curve 

The neutron derived porosity curve shown in Fig. 3 repre- 
that 


which had been caved with 7-in. pipe (Fig. 2). This curve 


sents data obtained from portion of the cored section 


may be applied to any other well in the area in which 7-in. 


casing is set and on which an absolute zero reference line 


was obtained 


rhis curve (Fig. 3). with its relatively small scattering, is 


typical of the reef limestones in this region and has heen 
found to be applicable over wide areas of West Texas. 

The same procedure was followed for the 614 in. open hole 
condition, The gamma ray and neutron curve deflections from 
zero were determined opposite each of the corresponding core 
porosities, and each deflection divided by its respective log 
thus The 


shown by the 


sensitivity converting to standard units. neutron 


standard units opposite clean limestone (as 


gamma ray curve) were then plotted as the asbscissa (linear 


cale) versus the per cent porosity as the ordinate (log scale) 


on semi-log graph paper and the average neutron derived 


porosity curve established by drawing a straight line through 
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these points. Fig. 4 shows the neutron derived porosity curve 
for 614 in. open hole. This curve may also be applied to any 
well in this area with 614 in. open hole on which a radio 


activity reference line was obtained. 


Comparison of Lansing - Kansas City Logs 
and Porosities 


Fig. 5 compares the logs and core porosity data for a well 
producing from the Lansing- Kansas City section, Stafford 
County, Kans. 

In addition, this example illustrates the application of a 
method of averaging the core porosities in order that the 
lengths of the core samples will agree more closely with the 
length of the sample affecting the neutron curve. The lengt! 
of the neutron curve sample is approximately two ft whil 
the core porosities were obtained from permeability plugs 
which were selected as representative of one-ft intervals of 
the core. This difference in length will result in a greater 
degree of scattering of the points used in determining the 
neutron derived porosity curve. The neutron curve subsurface 
instrument is not influenced to the same degree by each long 
itudinal unit of the interval sampled, as the middle portion 
of the sample (of approximately one ft of its total lengt! 
exerts the greatest influence while in the remainder of the 
sample (over the six in. at the ends) the influence of the 
sample on the neutron curve decreases to zero. The deter 
mination of the weighting factors would be too complicated 
to perform for a practical example and therefore a rough 
approximation was resorted too. This has been done by aver 
ing the porosities for three one-ft intervals of core so that 
the middle foot of core is weighted at full value while the 
foot below and the foot above is weighted by one-half: that 
is. the porosity of the middle portion is taken at full value 
one-half the value of the porosities of the adjacent intervals 
is added and their sum is divided by two. The result is shown 
in Fig. 5 in the column labeled “Per Cent Weighted Porosity. 
It will be noted that much better agreement exists betwee 
these values and the neutron curve than exists between th 
unweighted porosity values and the curve. 

The neutron derived porosity curve for the Lansing - Kan- 
City well is shown in Fig. 6. Here the neutron curve defle: 
tions are plotted in inches. 

To apply this neutron derived porosity curve to log 
other wells having the same physical bore hole conditions 
is necessary to apply a correction if the log was run or 
different This best be 


selecting two points on the base neutron derived porosit 


sensitivity. correction may made b 


curve, on even porosity lines, i.e., two per cent porosity and 
20 per cent porosity. The neutron curve deflection in inches is 
determined for each of these points and each is multiplied 
by the ratio of the sensitivity of the 
old log. in order to obtain the correct deflection to be plotted 


new log to that of the 


along the corresponding porosity line. These points are then 
connected by a straight line. giving the corrected neutror 
derived porosity curve. which may then be applied to the 


neutron curve being considered. 


Application to Logs for Which No 
Instrumental Zero Is Available 


The neutron derived porosity curve established from a giver 
log may be applied to another log having only a shale refer 
ence line. by determining the location of the shale refe 


rence 


line in relation to the absolute zero reference line on the base 


log. and measuring this distance in inches. which figure is 
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divided by the sensitivity of the base log to obtain the dis- 
tance from shale to zero in terms of standard units. The loca- 
tion of zero may be placed on the new log by multiplying the 
sensitivity of the log by the number of standard units obtained 
and positioning the zero this distance to the left of the shale 
reference line (Fig. 1). The neutron derived porosity curve 
determined from the base log may then be applied as if the 
absolute known. The accuracy of this method of 


zero were 


RADIOACTIVITY 
LOG 





+~ 


| | i | | 
INSTRUMENT ZERO 


FIG. 2— COMPARISON OF THE RADIOACTIVITY 
POROSITIES IN KELLY-SNYDER FIELD, TEXAS 


LOG AND CORE 
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FIG. 3 — NEUTRON DERIVED POROSITY CURVE FOR 7” CASED HOLE 
KELLY-SNYDER FIELD 


establishing absolute zero is dependent upon the constancy 


with which the average shale line is established from well to 


well and the accuracy of the response of the neutron curve 


APPLICATION OF QUANTITATIVE METHODS 
OF INTERPRETATION TO SAND FORMATIONS 


It has been known for time that the gamma ray curve 


saritie 
can be used to estimate the relative shale content of the sanc 
and for this purpose a linear relationship may be assumed to 
between the deflection of the curve and the 
the sand. The is affected 
solely by the gamma radiation emitted by naturally occurring 


the rocks 


radioactive 


exist gamma ray 


shale content of gamma ray curve 
certain of 
rocks the 


constant concentration. The 


radioactive elements distributed among 


composing the earth’s crust. In some 


elements are present in a very 


radioactive content of shale. for example. is often sufficiently 
constant in a given production zone that the percentage of 


shale at a given depth may be estimated by comparing the 
gamma ray curve deflection at the given depth with the gamma 
shale. The net 


responds to lithology in a 


deflection at a 
other 


ray curve point opposite pure 


tron curve, on the hand 


much more complicated manner. This curve reflects lithologic 
terms of the 


conditions in eflect: produced on the neutrons 


emitted from the source by the hydrogen content of the rock 
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the cross sections for the scattering and capturing of neutrons 
of the elements composing the rock, and the atomic density 
of these elements in the rock. 
In clean (silt and shale-free) limestone or dolomite forma- 
tions the effects of the hydrogen content and that of the cross 
section for neutron capture act in combination to shift the 
neutron curve to the left. The seattering cross section and 
and shift 


simple case where the neutron curve 


density effects also act in combination the curve 


to the right.” In the 


deflections opposite mixtures of dolomite or limestone and 


liquid hydrocarbons or water were studied, it was found that 
relation exists.” However, this simple 


a simple exponential 


relationship will not apply te most sandstones because of 
their more heterogeneous chemical composition and the result- 
the 
requirements for calculating a neutron derived porosity for 
the 


to the three factors previously 


ing more complex effect on the neutron curve. All of 


sandstones from neutron curve are not known, but in 


addition enumerated there 
may be factors which depend on the chemical composition 
of the 
ire still 


used to determine porosity in sandstones. 


matrix and cementing material. Because these factors 


unknown, the neutron curve has not. as yet, been 


However, it has been found that in certain cases the neutron 


and gamma ray curves can be used in combination to assist 


in estimating the productivity of a sand section: i.e., whether 


the section is too shalv to produce or whether it will produce 
liquid 


This use depends upon the reliability of the 





FIG. 4— NEUTRON DERIVED POROSITY CURVE FOR 6's” OPEN HOLE, 


KELLY-SNYDER FIELD 
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FIG. 5— COMPARISON OF THE RADIOACTIVITY LOG AND CORE 
POROSITIES IN STAFFORD COUNTY, KANSAS. 


gamma ray curve for estimating the amount of shale present 
and the application of this curve in a compensatory manne 
to the neutron curve. Shale contains large proportions of 
chemically bound water and hydrocarbons; and therefore the 
neutron curve is affected by shale in a similar manner to 
that produced by its constituent substances when in the free 
liquid state. The gamma ray curve may be used to calculate 
which has a compensatory effect for shale on the 
neutron curve deflections, and with its aid it is 
possible to estimate the productivity of shaly sands from the 
curve. This factor has been called the gamma ray 


a factor 
sometimes 


neutron 


shale factor. 


Calculation of the Gamma Ray Shale Factor (G.) 


The gamma ray shale factor is an empirically derived fun¢ 
tion of the per cent total gamma ray curve deflection. The 
total gamma ray curve deflection is the distance between the 
deflection values for clean, silt-free, shale-free sandstones and 
the deflection values for the pure shale. The per cent total 
gamma ray curve deflection is obtained by dividing observed 
ray curve deflection at a given depth by the total 
gamma ray curve deflection ,and multiplying by 100. G. is 
expressed as a logarithmic function of the total gamma ray 
curve deflection (D,) as follows: 

G. ] log [1 + log (D, + 1)] dae (1) 
G, varies from 1 to 1.477 as D, varies from 0 to 100. A plot 
of the values of G, versus the values of (D, + 1) are shown 


gamma 


in Fig. 7. 
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Application of Gamma Ray Shale Factor 
To Sand Formations 


{ relative method using an interrelation factor calculated 
from the per cent total neutron curve deflections and the 
gamma ray shale number has been developed, making it pos- 
sible to predict the productivity of a zone within certain lim- 
itations. This interrelation factor has been called the neutron 
productivity number (V,). Using a method based on the use 
of the neutron productivity number, it has been possible to 
determine the type of production to be obtained from the 
zone in question in many instances. 

The neutron productivity number (V,) is determined in the 
following manner: 

N,; = GAD 4+ 1) a 3 6 Oe Ree eee 
where: 

D, is the per cent total deflection for the neutron curve and 
G. is the gamma ray shale factor, or: 

ie (D 1) [1 + log[1 + log(D,4+1)]] . . (3) 

The neutron productivity number should be used only in a 
relative manner by comparing the productivity numbers cal- 
culated from the curves of a well whose productivity is to be 
estimated, with the values of productivity numbers for wells 
producing from the same section. To establish neutron pro- 
ductivity numbers for producing sections of a field, drill stem 
tests and production tests diagnostic for individual sections 
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FIG. 6 — NEUTRON DERIVED POROSITY CURVE FOR 4'2” CASED HOLE, 
STAFFORD COUNTY, KANSAS 
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FIG. 7 — RELATIONSHIP BETWEEN GAMMA RAY SHALE NUMBER AND 
GAMMA RAY CURVE DEFLECTION 


must be compared with the productivity numbers calculated 
from the logs. The neutron productivity numbers calculated 
from the log of a well whose productivity is unknown, but 
which penetrates the same sections can be compared to the 
neutron productivity numbers for the known 
estimate made of the probable production to be expected 


wells and an 


Calculation of the Neutron Productivity Number 

The variables D, and D, are determined from their 
First 
mined for each curve, which 


res pet 

average shale lines are deter 
are labeled “100” and “0” on the 
gamma ray curve and the neutron curve, respectively’ A typi 


tive curves as follows. 


cal clean section is then selected on the gamma ray curve 
and a line is drawn through these deflections parallel to the 
shale line previously drawn. This line is labeled “0” and the 
deflection between the line marked “O” and “100” is known 
as the “total gamma ray curve deflection.” (D,,). A 
drawn through the deflections of the neutron curve opposite 


producing liquid hydrocarbon or 


line Is 
a clean zone known to be 
water, and it also is drawn parallel to the shale line. This line 
is labeled “100° and the distance between this line and the 
shale line of the neutron curve is called the “total neutron 


(D,.) 
The per cent total curve detlections D, and D,, are obtained 


curve deflection,” 


by dividing the curve deflections opposite a particular depth 
and zone. measured from the appropriate zero line, by the 


196 


PETROLEUM TRANSACTIONS, AIME 


THE QUANTITATIVE APPLICATION OF RADIOACTIVITY LOGS 


total deflection in inches of that curve. (the distance between 

the zero and the maximum lines). This is illustrated in Fig. 8. 

lhe second step in calculating \V,, is to determine the values 

of D, and D, by measuring the deflection of the particular 

zone desired and comparing this value with the total deflec- 

tion for the gamma ray and neutron curves. respectively. Thus, 
D 


D D i arc Ne ae a den Sa 


where 
D is the 
curve deflection. and D 


the zero line for the 
is the distance for total 
Similarly for the neutron curve: 


distance measured from 
gamma ray 
deflection of the curve 
D.. 
A= a ae as ee ee ee ee et 

D 
rhe logarithms of the values for D, are then determined from 
any table of logarithms and the values substituted in Equa- 


tion (09). 

V, may be calculated more easily, however. by using the 
nomograph shown in Fig. 9, for here it is only necessary to 
lay a straight edge across the nomograph, aligning the values 
of (D, +1) and (D, + 1) on their respective indices and the 
value for ', will be found on the V 
intersection with the straight edge. 


index at its point of 


Response of the Neutron Productivity Number 


Vy, is an index which may indicate the per cent of shale, 
liquid or gas. The limits described in the following illustra- 
tion must not be taken as generally applicable, because it was 
found that these limits apply only in partcular areas. It is 
necessary that values of V, be determined and compared with 
drill tests, 
particular are 


stem production tests or perforating data in the 
i or field being considered before the range of 
values may be known. From such data the following table 
was compiled: 
Values of 
l to 75 
75 to 168 
, to 500 


Condition 

Shale 

Liquid Filled Porosity 
Gas Filled Porosity 


Specific Applications of the Neutron 
Productivity Number 


lo illustrate the application of the neutron productivity 
wells in the Freites District. Venezuela 


and the Buffalo Field of East Texas were selected. These wells 


number in sand areas, 
are thought to have penetrated a relatively uniform porosity 


throughout the sand bodies 


well in the 
2.323 ft was 
534 B/D. The 
maximum value reached by \V,, in this interval on the log 

The 7.445 - 7.455 ft 

tested but is believed by the operator to contain gas. The 
V,. value 
opinion, The interval 7.480 
believed by the 
for N 

The values 
ing the 


The first example. shown in Fig. 8. is of a 
District 7.295 - 


perforated and on test flowed at the rate of 


Freites Venezuela. The interval 
was 


approximately 150 interval was not 


of nearly 300 calculated from the log confirms this 
7.600 ft tested but is 
The value of 120 


was not 
operator to contain fluid 
tends to confirm this opinion. 

of \,, were calculated from the log by establish- 
reference lines in the manner described previously. 


and applying D_ and D, to the nomograph 
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Also shown in Fig. 8 is the radioactivity log of the Lon¢ 
Star Producing Co. Mary Jones No. 4, Buffalo Field. The 
main Woodbine sand from 5,695 ft to 5,724 ft apparently is 
a uniformly clean (silt and shale-free) sand. This zone wa 
tested through perforations in the 51 -in. casing from 5.704 
These 
perforations were squeezed and the casing reperforated from 
5.706 ft to 5.714 ft. A completion test of these perforation 
yielded 81.81 B/D of oil with little gas. 


ft to 5,712.5 ft and produced gas with very little oil 


The values for D, and D, throughout the entire section wer 
determined, and \,, was calculated from the nomograph (Fig 
9). The gas-oil contact is placed at 5,704 ft with a value of 
160 \,, (table above) representing the approximate separation 
value. 


The Lone Star Producing Co. W. D. Recknor No. 1, Buffalo 
Field, is another example of the application of the table. In 
this well the 5%5-in. casing was perforated from 5,741 ft to 
5,745 ft resulting in the production of oil with a low gas/oil 
ratio. It will be noted that the N, value at the depth 5,742 ft 
is 92 and at 5,740 ft the V, value is 260, placing the gas-oil 
contact between 5,740 ft and 5,742 ft. Although the section 
at 5.742 ft is indicated by the gamma ray curve as shaly, the 
value of 92 after application of the V, curve would indicate 
that it should be productive of liquid. The low values of \ 
in the interval from 5,728 to 5,736 ft are probably due to 
the invasion of large quantities of drilling fluid. The sands 
indicated from 5.754 ft to 5,758 ft are known to be water 
bearing and were used as the 100 per cent line on the neutron 
curve. The values of V, were computed by the use of the 
nomograph as previously described, Fig. 9. Had the NV, 
not been computed for this log, it would have been difficult 


curve 


to place the gas-oil contact within a two-ft zone due to the 
influence of shale at exactly the point of contact. 

The values established for this field are usable over a lim 
ited area and should be improved with each subsequent well 


WELL IN NIPA FIELD 
FREITES DISTRICT 
VENEZUELA 
GAMMA RAY NEUTRON 
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that is drilled and tested. The method described here is par- 
ticularly applicable to the shaly sands in which the influence 
of shale reduces the intensity of the neutron curve and pre- 
vents the visual identification from the neutron curve of gas- 
liquid interfaces which is possible in clean sands. 


SUMMARY AND CONCLUSIONS 


The use of a device for determining the position of zero 
radioactivity and a sensitivity selector switch of greater versa- 
tility has been discussed and their application to actual logs 
of the limestones of West Texas and Kansas has been shown 
by illustration. Using the new instrumental zero reference line, 
the application of quantitative methods to the radioactivity 
logs will provide a means of estimating porosity. This is done 
by comparing the intensity of the neutron curve deflections 
with the porosities from the core samples and plotting this 
relationship to establish a neutron derived porosity curve. 
This curve may then be applied to unknown wells of similar 
hole 
obtained. The data shown in this paper are applicable to 
wells having 7-in. cased holes or 614-in. open holes in the 
reef limestone of West Texas and to 41%-in. cased holes in the 
Lansing - Kansas City Kansas. 


bore conditions on which radioactivity logs have been 


zone in 

\ relative method of utilizing an interrelation factor calcu- 
lated and neutron curves called the 
neutron productivity number has been discussed. The neutron 
productivity calculated for given depth intervals 
within the zone of interest will aid the interpretation of radio- 


from the gamma ray 


number 


activity logs in determining gas-liquid contacts in shaly sands 
and sandy shales. The use of the neutron productivity number 
is recognized as a relative method of comparison and is lim- 
ited in applicability to those areas and reservoirs in which 


the neutron curve is definitive of gas in the reservoir. Greater 
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FIG. 9 — NOMOGRAPH FOR COMPUTING THE NEUTRON PRODUCTIVITY NUMBER 


versatility and a wider range of applicability may be provided 
by the use of new radiation zero devices and it is hoped their 
use will result in an extension of quantitative methods for 
sand areas. 


ACKNOWLEDGMENT 


The authors wish to thank the Lane-Wells Co. and Well 
Surveys, Inc., for permission to publish this paper. They also 
wish to express their appreciation to the oil companies whose 
cooperation and assistance made this paper possible. The 
authors wish to acknowledge the assistance of the members 
of the staffs of the Lane-Wells Co. and Well Surveys, Inc.. who 
have aided in the preparation of the manuscript. 


198 


PETROLEUM TRANSACTIONS, AIME 


REFERENCES 


. Green, W. G., and Fearon, R. E.: Geophysics, (July, 1940), 


2. Pontecorvo. B.: Oil and Gas Jour., (Sept. 11, 1941), 32. 

3. Russell, W. L.: AAPG Bul. (Sept.. 1941), 25, No. 9. 

t. Jackson, W., and Campbell, J. L. P.: Trans. AIME, 165, 241. 

5. Campbell, J. L. P.. and Winter, A. B.: Oil Weekly, (Jan. 6. 
1946). 

. Mercier, V. J.: Oil Weekly, (Oct. 14 and 21, 1946). 

. Fearon, R. E.: Nucleonics, (June, 1949). 4, No. 6, 30-92. 

. lbid. | April. 1949), 4, No. 4, 67-75. 

. Bush, R. E.. and Mardock, FE. S.: 
189, 19. 


AIME, (1950). 
x ke * 


Vol. 192, 1951 


Trans. 





THE SOUTHWEST ANTIOCH GIBSON SAND UNIT 


T. E. OCKERSHAUSER, THE GLOBE OIL AND REFINING CO., OKLAHOMA CITY, OKLA., MEMBER AIME 


INTRODUCTION 


The Southwest Antioch Oil Field located in T2&3N - R2&3W, 
Garvin County, Okla., was discovered in February, 1946, by 
The Globe Oil & Refining Co. and The Vickers Petroleum 
Co. at their Melinda Gibson No. 1. Oil production was en 
countered at a depth of 6,525 ft in 14 ft of lower Pennsyl 
vanian sandstone. 

The producing sand, locally called the 
Deese,” is a sand bar deposit strung along the west flank of 
the Pauls Valley uplift in a northwest-southeast direction. The 
oil reservoir is defined to the east by the edge of sand depo 
sition along the uplift and to the west by the gradation of 
sand into shale. Within the Antioch Unit, the sand dips from 
150 to 400 ft per mile southwesterly and has a pay thickness 
varying from 10 to 55 ft and averaging 28 ft. 

Several areas of Gibson production have been developed. 
including the Katie pool on the south, the Elmore Area of 
the Southwest Antioch Gibson Sand pool, the Southwest Anti 
och Gibson Sand pool proper, the Southwest Maysville pool 
and the New Hope pool on the north. Production extends over 
a belt 22 miles long by about two miles wide and is continu 
ous except into the Katie and New Hope pools. It appears that 
drilling will soon prove these intervening areas productive 
To March 1, 1950, approximately 18,000 acres have 
developed with 450 Gibson wells. It is likely this count will 
ultimately reach 32.000 acres with 800 wells. 


Gibson” or “Third 


been 


RESERVOIR 


The natural production of oil from the Gibson reservoir is 
primarily due to solution gas expansion though gravity also 
plays a part. There is no indication of water drive or evidence 
that an initial gas cap existed. 

Gibson oil has a stock tank gravity of 43° API and was 
saturated at the original reservoir pressure of 2,925 psi. It 
contained 1,340 cu ft of gas per bbl in solution and had a 
high shrinkage of 70 per cent. Connate water amounts to 15 
per cent. Reservoir temperature is 130°F. 

Porosity averages 15.5 per cent and permeability ranges 
from less than one md to over 1.000. Based on measurement 


Mauscript received in the office of the Petroleum Branch, Aug. 30, 195 
Paper presented at the Petroleum Branch Fall Meeting in New Orleans 
La., Oct.4-6,1950. 
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of 400 core plugs from 15 wells, 78 per cent of the perme- 
abilities fell within the 10 to 1,000 md bracket and averaged 
167. 

Primary oil recovery is estimated at 140 bbl per acre-ft or 
22 million bbl for the present 137-tract unit. Due to unitization 
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FIG. 1 — MAP SHOWING OlJL FIELDS IN VICINITY OF THE SOUTHWEST 
ANTIOCH GIBSON SAND UNIT. 
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and gas injection, it is expected that recovery will be increased 
to 260 bbl per acre-ft or to a total of 41 million bbl — about 
doubled. 


UNITIZATION 


It appeared to oil operators early in the development of the 
Southwest Antioch pool, that oil recovery could substantially 
be increased under a program of unit operation and gas 
injection. This belief was verified by a study completed in 
February, 1947. It was realized that to put such a program 
into effect, the full cooperation of royalty owners was neces- 
sary and desirable. To this end a meeting of royalty owners 
was called and operators helped organize the Southwest 
Antioch Royalty Owners Association. 


THE SOUTHWEST ANTIOCH GIBSON SAND UNIT 


This association played an important part in creating and 
perfecting the unit. Their technical experts actively partici- 
pated in operator's meetings and had access to all informa- 
tion. To keep its membership advised, there were several 
meetings of the association in the towns of Antioch and Mays- 
ville. Operator representatives attended a number of these 


gatherings to discuss unitization. 


On Jan. 6. 
heard and approved the petition creating the Southwest Anti- 


1948, the Oklahoma Corporation Commission 


och Gibson Sand Unit without objection from royalty owners 
or oil operators. Operation of the 67 wells comprising the 
original unit was assumed by The Globe Oil & Refining Co. 


as Unit Operator on April 1, 1948. 




































































FIG. 2 — STRUCTURE MAP, GIBSON SAND, SWAGSU 
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The Antioch Unit has been enlarged three times without 
incident and has grown in size from 67 to 137 wells. A fourth 
enlargement is new being considered which will increase this 
number to approximately 180 wells. 

Under the “Plan of Unitization” adopted by the Southwest 
Antioch Unit, only producing 40-acre tracts are taken in 
The original participation formula based half of tract interest 
on net pay and half on ability to produce up to 150 B/D of 
oil with a maximum of 375,000 cu ft of gas. The formula has 
since been modified to include a reservoir pressure factor 


Operation 


Operating provisions of the “Plan of Unitization” stipulate: 
“The oil produced from the Unit Area shall be produced 


T.P. 3057 


from those wells in the Unit Area from which the same can 
be obtained with the smallest loss or dissipation of reservoir 
energy reasonably possible under practical operating con- 
ditions as they may exist from time to time;” and 


“Gas wells and wells which produce oil with gas/oil 
ratios found to be excessive in relation to the gas/oil ratios 
of other wells producing oil from the Unit Area shall be 
shut in or the production therefrom restricted in such man- 
ner as to make the most effective utilization of the gas 
energy of the reservoir reasonably possible under practical 


operating conditions as they may exist from time to time.” 


Production operations from the beginning have been con- 
ducted in accordance with these principles. During the first 
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FIG. 3 — THICKNESS MAP, GIBSON SAND, SWAGSU. 
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month of operation, as a result of selective oil production, 
average gas/oil ratio for the Unit Area was lowered from 
3,000 to 2,200 cu ft per bbl or about a third more oil was 
produced for the same volume of gas. 


In order to determine the amount of oil to be produced from 
each well and determine the distribution of gas to input wells, 
gas/oil ratios are taken monthly and reservoir pressures are 
measured three to four times a year. 


Injection 

Gas from the Unit is processed and compressed for injec- 
tion by the plant of the Warren Petroleum Corp.. e¢ al. Injec- 
tion of 9 MMcf/D of gas started during August, 1948, with 
one 600-HP compressor. A second 600-HP machine, which was 


THE SOUTHWEST ANTIOCH GIBSON SAND UNIT 


placed in operation during October, 1948, boosted input to 15.5 
MMcf/D. Present capacity of 25 MMcf/D was reached in 
June, 1949, with the addition of a third 800-HP compressor. 

Gas is being returned to the Gibson sand through seven 
injection wells situated at high structural positions near the 
sand edge on the east side of the pool. Injectivity rates at 
individual wells range from 2 MMcf/D to in excess of 5.5 
MMcf/D at 2,400 psi pressure. During 21 months operation, 
a total of over 10 billion cu ft of gas has been put back into 
the reservoir 


The effects of gas injection were first noticed during the 
second month of injection when a material lessening in the 
rate of bottom hole pressure decline occurred near input wells; 
some wells even showed small pressure increases. At one pro- 
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FIG. 4 — BHP CHANGE MAP, SEPTEMBER, 1948, THROUGH FEBRUARY, 1950, SWAGSU. 
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ducer, during the third month of injection, oil rate tripled 
and gas/oil ratio fell to one-third. During the fourth month 
at another well, oil gravity increased from 40° to 60° API. 
This change in liquid character was caused by injection of 
gas at a diagonally offsetting well which formerly produced 
high gravity oil. 


The advance of the gas front away from injection wells has 
been fairly uniform, conforming closely with structure. There 
has been no tendency for gas to channel or break through. 
From a comparison of pressure drops at individual wells with 
the average drop for the Unit since the beginning of injection, 
it appears gas has influenced 45 wells, an area of about 2,000 
acres, and has been felt as far away as 6,500 ft from input 
wells. 


PERFORMANCE 


The 137 wells comprising the present Unit have produced 
a total of 10.8 million bbl of oil to March, 1950. Of this 
amount 9 million bbl or about 85 per cent has come from the 
80 wells in the first enlarged Unit. Inasmuch as a longer 
history is available on these wells and they have recovered 
the major portion of the oil, their behavior is a reliable indi- 
cation of reservoir performance. 


At the time of unitization, April 1, 1948, gas/oil ratio was 
2,900 cu ft per bbl and was increasing rapidly. Had conven- 
tional operation continued, gas/oil ratio by now would be in 
excess of 10,000. Actually it is 4,800. It appears that a mini- 
mum of twice as much oil is now being recovered per cu ft 
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FIG. 5 — TIME-POSITION 10,000 GOR CONTOUR, VICINITY GAS INJECTION, SWAGSU. 








Vol. 192, 1951 


PETROLEUM TRANSACTIONS, AIME 





T.P. 3057 


saso 














POUNDS PER SQUARE INCH aT 





RLEERVO PRESSURE 














SESS SSS 
CUMULATIVE PIPE 
e 


[UNE RUNS (MILUONS OF BARRELS) 


FIG. 6 — RESERVOIR PERFORMANCE, 80 TRACT UNIT, SWAGSU 


than would have been recovered had conventional 


operation continued. 
\ total of 5.2 million bbl of oil and 18.1 billion cu ft of 
gas have been produced under Unit operation and 9.2 billion 


ot gas 


cu ft of gas have been returned to the reservoir. Net gas/ oil 
ratio has averaged 1,725 cu ft per bbl. This compares to an 
average ratio of 1,645 cu ft per bbl before unitization, during 
which period 3.9 million bbl of oil and 6.4 billion cu ft of 
gas were produced. Considering only the interval of gas injec- 
tion, 4.2 million bbl of oil have been recovered at a net ratio 
of 1.565 cu ft per bbl. 

The effect of gas injection is readily apparent from an ex- 
omination of reservoir pressure decline which shows a decided 
reduction in rate during December, 1948, five months after 
gas injection started. Prior to this time 7.370 bbl of oil were 
produced per pound pressure drop. Since this date, each 
pound drop in pressure has accounted for the recovery of 
17.790 bbl or approximately 140 per cent more oil. 

The public, much to the harm of unitization, frequently is 
of the impressian that reduced production rate inevitably ac- 
companies unit operation. Comparison with the production rate 
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of competitively operated wells in the Elmore Area of the pool, 
shows that such a belief is unjustified, at least in the case of 
the Antioch Unit. 

Both groups of wells produce from the same source of sup- 
ply, wells are very nearly the same age, and the reservoir is 
approximately in the same state of depletion. The Unit has 
produced the top allowable rate set by the Corporation Com- 
mission; whereas, the average well in the Elmore Area has 
produced substantially under the maximum allowable. Over 
a period of 17 months since gas injection started, Unit wells 
produced an average of 93.5 B/D of oil per well compared to 
73.8 B/D for wells in the Elmore area. The difference of 19.7 
B/D amounts to the production of 7,200 bbl of additional oil 
in a year’s time from each well in the Unit. To royalty owners 
this means an increase in annual income of $2,400.00 per well. 
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FIG. 7 — COMPARATIVE APPROACH OF SWAGSU AND ELMORE AREA 
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STEADY FLOW OF TWO-PHASE SINGLE-COMPONENT 
FLUIDS THROUGH POROUS MEDIA 


FRANK G. MILLER,* U. S. BUREAU OF MINES, SAN FRANCISCO, CALIF., MEMBER AIME 


ABSTRACT 


This report presents developments of fundamental equations 
for describing the flow and thermodynamic behavior of two- 
phase single-component fluids moving under steady conditions 
through porous media. Many of the theoretical considerations 
upon which these equations are premised have received little 
or no attention in oil-reservoir fluid-flow research. The signifi 
cance of the underlying flow theory in oil-producing operations 
is indicated. 

In particular, the theoretical analysis pertains to the steady, 
adiabatic, macroscopically linear, two-phase flow of a single 
component fluid through a horizontal column of porous me 
dium. It is considered that the test fluid enters the upstream 
end of the column while entirely in the liquid state, moves 
downstream an appreciable distance, begins to vaporize, and 
then moves through the remainder of the column as a gas- 
liquid mixture. The problem posed is to find the total weight 
rate of flow and the pressure distribution along the column 
for a given inlet pressure and temperature, a given exit pres 
sure or temperature and given characteristics of the test fluid 
and porous medium. 

In developing the theory, gas-liquid interfacial phenomena 
are treated. phase equilibrium is assumed and previous the 
oretical work of other investigators of the problem is modified 

Laboratory experiments performed with specially designed 
apparatus. in which propane is used as the test fluid, substan 
tiate the theory. The apparatus, materials and experimental 
procedure are described. Comparative experimental and the 
oretical results are presented and discussed. 

It is believed that the research findings contributed in this 
paper should not only lead to a better understanding of oil 
reservoir behavior, but also should be suggestive in regard to 
future research in this field of study. 


INTRODUCTION 


In recent years much time and effort has been 
in both theoretical and experimental studies of the static and 
dynamic behavior of oil-reservoir fluids in porous rocks. Al 
though lack of sufficient basic oil-field data, principally con 
cerning the properties and characteristics of reservoir rocks 
and fluids, largely precludes quantitative application of re 
search results to oil-field problems, qualitative application has 


In effect. 


consumed 


become common practice. oil-reservoir engineering 


Manuscript received in the office of the Petroleum Branch Sept 
1950. Paper presented at the West Coast Fall Meeting of the Petroleum 
Branch, AIME, in Los Angeles, Calif., Oct. 12-13, 1950. 
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research is serving as a firm foundation for oil-field develop- 
ment and production practices leading to increased economic 
f petroleum. This province of research, however, 
still poses many perplexing problems. 


recover ies 


The thermodynamic behavior of two-phase fluids moving 
through porous media constitutes one facet of reservoir-fluid- 
flow research that has not received the attention it deserves. 
This report embodies a theoretical discussion of this subject 
and a description of a series of related laboratory experiments. 
The significance of the problem to oil field operations is indi- 
cated but in particular the report centers around a theory and 
method for analyzing the steady, macroscopically linear, two- 
phase flow of a fluid (a single molecular species) through a 
horizontal column of porous medium. 

For simplicity in showing how the thermodynamic behavior 
of two-phase fluids moving through porous media affects oil- 
performance problems, attention is focused tempo- 
rarily on a particular well producing petroleum from an 
idealized water-free solution-gas drive reservoir, the reservoir 
rock being a horizontal, thin, fairly homogeneous sandstone of 
areal extent confined between two impermeable strata. 


reservoir 


large 

The flowing hydrocarbon fluid is considered to exist entirely 
as a liquid at points in the reservoir remote from the well; 
decline in fluid pressure in the direction of the 
vaporization of the hydrocarbon to begin at a 
radial from the well. Upstream from r the fluid 
moves entirely as a liquid and downstream from r it moves 
gas or as a gas-liquid mixture depending 
hydrocarbon and on the thermo- 


however, the 
well causes 

distance r 
either entirely as a 
on the 
dynamic process it follows during flow. 


properties of the 


The distance r would be variable under transient flow con- 
ditions, but for purposes of analysis the flow is considered to 
be oteady at the particular instant of observation during the 
of the well of interest. If the flow were isothermal 
hydrocarbon a pure substance, the fluid would be 
entirely gaseous downstream from r. Thus, this isothermal flow 
process for a pure substance would require that the heat of 
vaporization be supplied at r. over zero length of porous me- 
dium. at the precise rate necessary to maintain the constant 
temperature. This means that the solid matrix of the porous 
medium (reservoir rock) and the surroundings (impermeable 
strata confining the reservoir rock) would have to serve as 
heat requirements would be 
isothermal flow of a multicom- 
ponent hydrocarbon as bubble and dew points at the same 
temperature correspond to different pressures. In this instance 
isothermal conditions would be sustained without complete 
vaporization of the fluid over zero length of porous medium. 


Nevertheless, 


flow ing life 
and the 


infinite sources. Heat-transfer 


somewhat less severe for the 


as the flow is in the direction of decreasing 
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pressure, isothermal flow of either a single-component or multi- 
component hydrocarbon represents a maximum limiting condi- 
tion as no other thermodynamic, decreasing-pressure, vapori- 
zation process would require as great a rate of heat transfer 
to the fluid. In contrast, a minimum limiting condition is rep- 
resented by a thermally-isolated flow system. defined as one 
in which no heat is transferred to the flowing fluid from the 
solid matrix of the porous medium or from the impermeable 
strata confining the system. 

As pointed out by O’Brien and Putnam’ the flow conditions 
in an oil reservoir apparently are intermediate between the 
foregoing maximum and minimum limiting conditions. Those 
authors indicate that a petroleum-bearing stratum could be 
represented by a fluid-filled porous medium, confined above 
and below by constant-temperature impermeable beds from 
which heat is transferred along various paths, solid and fluid, 
through the porous media. Flow problems relating to this 
system, however, depend on complete knowledge of the heat- 
transfer characteristics of the system and such information is 
not presently available. In this connection O’Brien and Putnam 
indicate further that isothermal and adiabatic flow may be 
very nearly the same for the flow of a crude oil and natural 
gas mixture in an oil reservoir. They cite an example pertain- 
ing to the adiabatic flow of such a mixture and show that the 
fluid temperature drops only a few degrees for a drop in pres- 
sure from 3,000 to 400 psia. 

The laboratory study described in this report is based on 
the foregoing considerations. In an effort to start with the 
simplest appropriate two-phase flow system it was deemed 
advisable to investigate the two-phase flow behavior of a 
single-component fluid. This restriction rules out isothermal 
flow as a single component will not flow isothermally in two 
phases. Accordingly, adiabatic conditions were adopted. 

In studying the adiabatic system it was found necessary to 
modify and extend the O’Brien-Putnam flow theory. The 
modified theory, concerning fundamental relations between 
rates of flow and pressure differentials, is substantiated by 
the experiments. 


SCOPE OF PROBLEM 


The investigation described in this paper is intended to con 
tribute to the solution of the general problem of determining 
the flow behavior of oil-reservoir fluids in subsurface porous 
rocks. It is well known that this general problem represents 
a large field of widely diversified research in the petroleum 
producing industry. The scope of the present work is con 
fined to the development and test of equations that describe 
the flow and thermodynamic characteristics of a single-compo- 
nent two-phase fluid moving through a porous medium under 
conditions of adiabatic flow. 

The fluid-transmitting capacity 
porous medium have no direct bearing on the principles of 


and the geometry of the 
the analysis. For convenience. therefore. the test flow system 
was made by packing a suitable conduit with a clean, previ 
ously sized unconsolidated sand. As the bulk cross-sectional 
area of the resultant sand column was sensibly constant, the 
flow was macroscopically linear. Propane was selected as the 
test fluid. 
available for this substance. 
was selected in preference to other fluids because of the inter 


Adequate thermodynamic and viscosity data®’ are 
Moreover, a hydrocarbon fluid 
est of the oil-producing industry in the behavior of hydrocar- 
bons. To simplify the analysis the flow was specified as steady 

In’ conformance with the preceding discussion the flow sys 
tem was designed to prevent transfer of heat between the 


flowing fluid and its surroundings (the sand grains and the 
1References given at end of pape 
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walls of the flow conduit). Furthermore, in the experiments, 
the propane entered the column of sand as a liquid, began to 
vaporize an appreciable distance from the inlet and discharged 
from the column as a gas-liquid mixture. 

To avoid the introduction of unnecessary complications 
regarding the meaning of pressures measured in the two-phase 
flow portion of the system, the column of sand was installed 
in a horizontal position. 

Five flow experiments were performed. Total 
liquid) weight rates of flow were measured and pressure and 


(gas plus 


temperature distribution data were collected. 

The ensuing theory, substantiated by the propane-flow ex- 
periments conducted with the described flow system, discloses 
that gas and liquid weight rates of flow as a function of dis- 
tance along the sand column, and pressure and temperature 
distributions, calculated if the inlet pressure and 
temperature and the outlet pressure or temperature are given. 
In addition, pressure and temperature distributions and gas 
and liquid weight-rate-of-flow distributions, may be calculated 
if the constant total weight rate of flow and the inlet pressure 


may be 


and temperature, are given. 

In developing the theory, capillary phenomena are treated. 
Based on the work of Aboul-Seoud,’ it is assumed that phase 
equilibrium is attained at all cross-sections of the column of 
sand perpendic ular to the direction of flow. 


THEORY 
Nomenclature 
The symbols used in mathematical analyses in this report 
follows: 
x = distance along column of porous medium 


are as 


value of x at which vaporization begins 
subscript denoting liquid 
subse ript denoting gas 
weight rate of flow (with no subscript, W = total 
+ W,) 
- enthalpy per unit weight of fluid 
kinetic 
= fluid pressure 
= fluid pressure at upstream end of column of po- 


weight rate = WV, 
energy per unit weight of fluid 
rous medium 


fluid 
porous medium 


pressure at downstream end of column of 
equilibrium, plane-surface, vapor pressure 
absolute temperature 

weight of gas per unit weight of gas-liquid mix- 
transverse section of column in 


ture passing a 


unit time 


weight of gas per unit weight of gas-liquid mixture 
instantaneously at a transverse section of column 
(with no subseript, Q = 
+ Q,) 


volumetric rate of flow 
total volumetric rate = Qy 
specific volume 


bulk 


medium 


cross-sectional area of column of porous 
fractional porosity 
wetting-phase content (fraction of pore space oc- 
cupied by liquid) 
= time variable 
average miscroscopic velocity 
homogeneous-fluid permeability 
permeability to liquid (with gas present) 
permeability to gas (with liquid present) 
absolute viscosity 
Fy (p) ky/h 
fk {p) a h 
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ky /k 
pressure function = : — + 
My Vi 


. ke/k 
F(P) —- 
MeV 
P. = capillary pressure 

R = universal gas constant 


M = molecular weight of gas 


Development of Flow Equations 


For steady two-phase (gas-liquid mixture) flow, continuity 
principles immediately lead to 
dW dW, 
~- — = 0 a eee 
dx dx 
Vo+ = F m= ecomtemt.....« .. 
The general energy equation, which is premised on condi 
tions of steady flow, provides means for determining the ther 
modynamic process undergone by the flowing fluid. Regarding 
that equation, it is evident that the described flow system 
involves no “shaft work” and that the horizontal position of 
the column implies no changes in potential energy. Therefore 
ae a es a oS oe et 
Equation (3) is not limited to single-phase flow; it applies 
equally well to two-phase flow. In either instance, however, 
H and U are to be based on a unit weight of fluid flowing. 
If it is assumed that Darcy’s law is valid for the system 


under discussion, viscous forces would predominate and the 


heat energy generated by friction would be very large com 
pared to kinetic energy changes. Thus, for any two transverse 
sections in the system, a finite distance apart, 

PS Ee ee ee ee eee |, 
hence, 

H=constant. . - « 


which means that the fluid follows an isenthalpic process. As 
a consequence the fluid will flow in two phases (gas and lig 
uid) if the (H, T) or (H, P) conditions along the column 
fall within the two-phase region on an ordinary thermodynamic 
chart for the test fluid (Fig. 1). 

Based on Equations (2) and (5) the enthalpy of the fluid 
passing any transverse section in the column in unit time is 


given by 
W.H, + W.H, = WH=constant .... (6) 
Suppose that the fluid enters the column (at x = () entirely 
as a liquid but that the flowing-fluid pressure has become re- 
duced to its saturation intensity at x». For an isenthalpic prox 
ess the fluid must exist downstream from xo as a gas-liquid 
mixture. Under these conditions the flowing fluid would follow 
a path as represented by path ABC in Fig. 1. 





ENTHALPY 














L 
TEMPERATURE 


FIG. 1— SCHEMATIC ENTHALPY-TEMPERATURE CHART FOR PROPANE 
(AFTER O'BRIEN AND PUTNAM). 
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In unit time a weight W of liquid enters the column, and 
during the same unit of time a weight W,, of liquid and VW, 
of gas departs from any arbitrarily selected transverse section 
at x. Thus, for the length of column bounded by zero and x, 
the fraction (f) that vaporizes, of an original unit weight of 
entering liquid, may be expressed by 

N 
Substituting in Equation (2 


j= 


hence, 


Equation (6) becomes 


(1-f)H, + fHe =H (10) 


f= ive. oi aes tk eee ta, Sa 


In the thermodynamics of single-component static systems, 
the righthand member of Equation (11) is a well-known 
expression for “quality,” the weight of gas present in a unit 
weight of a gas-liquid mixture. For a static system H would 
be the enthalpy per unit weight of mixture present, and H, 
and H,, the respective enthalpies per unit weight of gas and 
liquid present. For the foregoing dynamic system (see Equa- 
tion 6), it has been shown that an isenthalpic process is one 
which requires the enthalpy of the fluid passing any transverse 
section in unit time to be the same for all transverse sections 
of the column. Consequently, for a dynamic system enthalpies 
are to be expressed on a per-unit-weight-of-fluid-flowing basis. 
As charts like Fig. 1 are prepared from data collected with 
static systems, it is desirable to show how they may be applied 
to dynamic systems of the kind under discussion. 

In this respect. suppose that it is possible to separate and 
collect the gas and liquid discharged in unit time from the 
transverse section at x. If the collected gas is kept entirely 
in the gaseous state at the pressure and temperature at x, it 
will be saturated and occupy a volume Qo, numerically equal 
to the volumetric rate of gas flow at x. Likewise, if the col- 
lected liquid is kept entirely in the liquid state at the pressure 
and temperature at x, it will be saturated and occupy a volume 
Q,, numerically equal to the volumetric rate of liquid flow at x. 
Furthermore, if these collected and separated fluids are recom- 
bined to the total volume (Q, + Q,.) without change in their 
common pressure or temperature, the weights of gas and liquid 
in the resultant mixture will be equivalent numerically, to the 
respective weight rates of flow of the two phases at x, namely 
W., and W,. The quality of the resultant mixture, therefore, 
will be W/W. 


Accordingly, the quality of the mixture passing x in unit 
time is W/W. which is f (Equation 7). Hence, the lines of 
constant quality on Fig. 1 are lines along which the fraction 
vaporized (f) is constant. 

(gain considering the foregoing hypothetically collected and 
recombined fluid, produced from transverse section x in unit 
time. the specific volumes of the collected gas and liquid, 
. are equal to Q./W. and Q,/W,.; thus. 
the specific volume of the resultant mixture as a whole is 
given by 


respectively ¥, and 7 


Q 
;, Regehr tl a 
ma ae W —_ 


It follows that 7 is the specific volume of the mixture passing 
x in unit time. 
It is to be emphasized that the weight fraction of gas in 
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the mixture instantaneously at x is quite different than the 
weight fraction of gas in the mixture passing x in unit time. 
In this respect the weight of gas instantaneously in an ele- 
mental length 4x of the column, upstream but adjacent to x, 
is (1/%q) [A.4x.€. (1—p)] and the weight of liquid is 
(1/¥,) (4. dx. €.) so that the weight fraction f, of gas 


in the element is 


1 


+(*)F 


To find a corresponding expression for the weight fraction 
of gas in the mixture passing x in unit time, attention is to 
cused on average microscopic velocities of the two phases. In 


j= (13) 


this connection, the liquid is assumed to be the wetting phase 
and the gas the nonwetting phase. Part of the gas may be 
stationary in the form of discrete bubbles enlodged in pores 
surrounded by liquid. Part of the liquid may be stationary 
owing to the blocking of its passage through gas-locked pore 
channels. Neglecting pressure differentials, at x, across curved 
gas-liquid interfaces, the enthalpies of the flowing gas and 
stationary gas are the same because they are subjected to the 
same pressure and temperature. Similar reasoning applies to 
the liquid. 

Let the average particle of liquid at the upstream end of 
the foregoing elemental length of column reach station x in 
time 6f, and the average particle of gas reach station x in 
time dt’. Then, as 6x may be made arbitrarily small, the aver- 
age microscopic velocity of the liquid in the macroscopic 
direction of flow, at x, is (4x/ét) likewise. the corre 
sponding microscopic velocity of the gas is (4x /5t’) 


= 3 
= Ma. 
Based on continuity principles 


(l—p)uy 


hence, 


We 
We + 


] 
R Vo p uy; 
1 + 
7 l-p Fu 
Comparison of Equations (13) and (16) indicates that / 
and u, are equal. The relationship 


and f are equal only if u 
between these microscopic velocities can be investigated by 
application of Darcy's law to the two phases independently* 
and by recourse to Equations (14) and (15): 


oot aes 1k dP = 
I . , (17) 
y dx 


dP 


° (18) 
dx 


W 


from which 


(19) 


. (20) 
F.,(p) u 

Knowledge of fluid viscosities and of the characteristics of 
relative permeability curves (for example. see References 5 
6 and 7). shows conclusively that it would be wholly fortuitous 
for the righthand member of Equation (20) to equal unity 
hence, u, and u, cannot be assumed equal. Consequently, it 
would be invalid to assume equality of f and f as was done 


*Wrebol and Botset e the first z to apply Darcy 


in two-phase flow 
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by O’Brien and Putnam (X in their nomenclature is the same 


as f in this report). 
Applying Darcy’s law, Equation (2) becomes 


ky /k keg /k dP 
It 1k 4 : 
y dx 


Integrating over a length of column between the inlet end of 


(21) 


My Y Maa 


the column and an arbitrarily selected transverse section at x: 
2 


I 
Wx th f 
P 
in which F(P), the pressure function, is equivalent to the term 
in brackets in Equation (21), and P, the fluid pressure at the 
inlet 


F(P)dP piace gs edie gy eee 


Equation (22) to determine the theoretical 
from the system. In this instance let x be 


the total length of the flow column. causing P to be the outlet 


may be used 


discharge rate, W. 


pressure. If small increments of pressure are taken along path 
{BC in Fig. 1, then Equation (11) may be used to determine 
the distribution of f along the column. O’Brien and Putnam 
advanced this method of calculation, but in using f as though 
setting / 
thereby 


it were the same as f they applied Equation (13) 
equal to the righthand member of that equation 
calculating p erroneously. They then used the incorrect values 
conjunction with relative-permeability curves to ob- 
k and k,,/k. Thus, their theory requires modification. 
In this respect, it is only necessary to use Equation (9) for 


k and k,,/k. From Equa- 


of pu 
tain k 


determinations of distributions of k 
(17) and (18): 

W A yu 

Kt ko 
Substituting for (B n 


tions 
(23) 
) in Equation (9) and rearranging 
; iV 1-f 
aera Gs, 


can be used to determine the k,/kg 


terms 
(24) 


Equation (24) ratio 
distribution along the column. 
data show k,/k. ke/k and kg/k,, as 
functions of with the k,./k, distribution known, it 
is possible to plot F/P) against P. The area under the result- 
vetween the limits of P 
of the integral in Equation (22). It 


of relative permeabilities 
Relative-permeability 


hence 


and P represents the value 


follows that 


ing curve 

may be 
calculated. 
With the 

» 


lished. Equation (22 


theoretical total weight rate of discharge estab- 
) may be applied further to determine the 
theoretical pressure distribution. To accomplish this calcula- 
tion, successive values of P (less than the inlet pressure and 
greater than the outlet pressure) are selected. In each instance 
(22) is evaluated graphically and 


the integral of Equation 


corresponding values «f x determined 


Gas-Liquid Interfacial Curvature 

Charts like Fig. 1 are based not only on equilibrium thermo- 
dynamics, but are. constructed from data collected with instru- 
ments that cause gas-liquid interfaces to be. sensibly. plane 
-urfaces. In the foregoing analysis, therefore, it is assumed 
that equilibrium thermo-dynamics applies and that the curved 
surfaces of gas-liquid phase separation in the porous medium 
may be treated, thermodynamically, as plane surfaces. 

According to the kinetic 


change of molecules between the phases of a single-component 


theory there is a continuous ex- 


gas-liquid mixture, phase equilibria being attained when the 


A boul- 


Seoud’ has made a theoretical proof, substantiated by experi- 


rate of vaporization and rate of condensation are equal. 
ment, indicating that phase equilibrium is attained almost 


instantaneously in a single-component two-phase flow system 


of the kind under discussion. 
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It is desirable, however, to indicate the effects of gas-liquid 
interfacial curvature and disclose a means for determining 
whether capillary pressures are negligible. By definition 

PPE, OLS ie or ee ae en's 
hence, the gradient of the capillary pressure in a direction 2 
is given by 

dP. dP, dP, 

—— = -— ete Hen See 

dx dx dx 

It is pointed out that Equations (25) and (26) have been 
derived for static systems only. As a first approximation, there 
fore, suppose that Equation (26) applies to the gas-liquid 
mixture instantaneously at an arbitrarily selected transverse 
section in the two-phase flow portion of the test column of 
porous medium. Then, the pressure gradients in Equations 
(17) and (18) would differ by (dP./dx). 

Considering that the liquid and gas respectively are the 
wetting and nonwetting phases, the Kelvin-Helmholtz equa- 
tion,”” pertaining to static equilibrium at curved surfaces of 
gas-liquid phase separation, may be written 

RT P,’ i 

P. = Pg-P.’ + ln ae te ae ae 4, 

Vi, Po 
(Assuming that this equation applies at any arbitrarily selected 
transverse section in the two-phase flow portion of the column 
of porous medium, and substituting for P. from Equation (25). 
results in 


(25) 


RT P,’ 
Pi = Pe ln 

Vy, Po 
{s the liquid is the wetting phase P,,’ > Pg; hence, Py < P.' 
Therefore, if liquid pressures and temperatures are measured 
along the column of porous medium in the direction of flow 
and in the two-phase flow portion of the system, gas-liquid 
interfacial curvature would be manifest by liquid temperatures 
greater than the normal (plane-surface) saturation tempera 
tures corresponding to the measured pressures. 

The same line of reasoning indicates that gas-liquid inter- 
facial curvature is unimportant, thermodynamically. if the 
corresponding values of the liquid pressures and temperatures 
are plotted on a pressure-temperature diagram and found to 
coincide, within the limits of experimental accuracy, with 
the normal {plane-surface) equilibrium vapor-pressure curve 
Under these conditions. P,, = P,’. and from Equation (28) it 
follows that P,’ = P.: thus P, = 0 and (dP./dx) = 0. 

Equation (22) was developed without consideration of a 
capillary-pressure gradient distribution (Equation 26) along 
the two-phase fluid portion of the flow column. Experimental 
verification of the flow theory — applying Equation (22) 
will depend, therefore. on whether capillary pressure gradient- 


are negiligible. 


ZTHERMOCOUPLE LEAD WIRES 





eT [ a oe = 
oie ae * 








FIG. 2— SCHEMATIC DIAGRAM OF APPARATUS USED IN PROPANE 
FLOW EXPERIMENTS. 
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In conformance with Muskat’s” reasoning pertaining to the 
characteristics of “capillary pressure-desaturation” curves, the 
existence of the “equilibrium gas saturation” will preclude 
any discontinuity in capillary pressure at the transverse sec- 
tion in the flow column where vaporization begins. Recognizing. 
moreover, that vaporization is a continuous process throughout 
the two-phase fluid portion of the system, there would be no 
discontinuous accumulation of the wetting phase at the outlet 
end of the column. Consequently, if, also, there are no strati- 
fication discontinuities in the porous medium comprising the 
column, there will be no discontinuities in the distribution of 
capillary pressure along the column. Thus, negligible capillary 
pressure may be taken as tantamount to negligible capillary- 
pressure gradients. 

For a single-component liquid in the vicinity of the satu- 
rated liquid line (path AB, Fig. 1), (@7/@P) « usually is quite 
small. For this reason the isenthalpic process followed by the 
flowing liquid in the all-liquid flow portion of the system will 
very nearly coincide with an isothermal process. If, therefore, 
a normal vapor-pressure chart is plotted and pressure and 
temperature distribution data, collected from the experimental 
flow system, are superimposed upon the chart, the all-liquid 
flow portion of the system should be indicated by a series of 
points forming a straight line parallel to the pressure axis, the 
pressures corresponding to these points being greater than the 
saturation pressure corresponding to their common tempera- 
ture. Points based on data collected in the two-phase fluid 
portion of the system will fall, within the limits of experi- 
mental accuracy, on the vapor-pressure curve, if capillary 
pressures are negligible. This procedure for plotting data was 
followed with respect to the ensuing propane-flow experiments. 
The results indicate isothermal flow in the all-liquid flow 
portion of the column of sand and that capillary pressures 
are small. 

It is of further significance to note that the foregoing straight 
line through pressure-temperature points representing all- 
liquid flow, may be extrapolated to the vapor-pressure curve 
to find the pressure and temperature at which vaporization 
begins. These data together with the experimentally-determined 
pressure and temperature distribution curves will serve to 
indicate the transverse section at which vaporization begins. 


EXPERIMENTAL EQUIPMENT AND 
TECHNIQUE 


Apparatus 


The apparatus used in making the experimental study of 
the adiabatic, two-phase. steady flow of propane through a 
horizontal column of sand is illustrated schematically in Fig. 2. 
rhe flow tube was about five ft long and was constructed from 
six lengths of nominal two-in. diameter “extra-streng” steel 
pipe, connected together in a line by forged-steel flanges. 

The flange-type units at the extreme ends of the sand-packed 
conduit were made by drilling. tapping and otherwise prepar- 
ing blind flanges to provide the confined column of sand with 
an entry and exit for the test fluid. Blind flanges also were 
used to construct the flange-type units to which manometer 
lines and thermocouple lead wires are connected (Fig. 2). 
Each of these special units was equipped with a built-in ring- 
type piezometer, similar to the kind made by Wyckoff and 
Botset.” and was machined to accommodate a fluid-tight, elec- 
trieally insulated coupling device. which permitted insertion 
of the thermocouple (copper-constantan) in the 
system and in the sand, without the necessity for conventional 


thermocouple wells. 
The flow tube composed of the pipe sections and flange units 
with a layer of heat-insulating material (mag- 
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nesite) about one-half in. thick (this coating is not shown in 
the sketch). Furthermore, the tube as a whole was enveloped, 
as shown, by a circular, compartmented sheet-metal air duct. 
Hot and cold air lines leading to each compartment provided 
the means for fixing the temperature distribution in the air 
along the outside of the flow tube to correspond to the tem- 
perature distribution of the flowing propane inside the tube. 
In this respect, the thermocouple circuit (not detailed in the 
sketch) that was used, was arranged so that temperature dif- 
ferences between the outside and inside of the flow tube could 
be detected. 

The separating chambers were made from nominal 1%%-in. 
diameter steel pipe and were about two in. long. They were an 
integral part of the manometer lead lines and were mounted 
on the sand-packed flow tube, inside the sheet-metal duct. 
During a two-phase flow experiment and while pressures were 
being measured, the propane in each chamber along the two- 
phase flow portion of the column of sand existed as a saturated 
liquid in connection with and at the temperature of the fluid 
inside the main conduit. Liquid water separated the propane 
in the chambers from the mercury in the manometers. 

The pressure gauge at the inlet end of the system was an 
accurately-calibrated test gauge having a maximum capacity 
of 300 psi and a least count of one psi. 

Each of the manometers was a glass U-tube assembly long 
enough to indicate a maximum deflection of 48 in. in the mer- 
cury levels in the two limbs and strong enough to withstand 
a total fluid pressure of at least 265 psia. All the manometer 
lead lines were copper tubing. equipped with brass needle 
valves and Imperial fittings. 

The first manometer on the left (Fig. 2) indicated the pres- 


sure drop across the inlet to the flow system. The entering 
liquid first passed through a Tyler 200-mesh copper retaining 
screen, against which the sand was packed, and then diverged 
to fill the sand column. Similarly, the last manometer on the 
right indicated the pressure drop across the outlet. which also 


was provided with a retaining screen. and which caused con- 
vergence of the flow lines. The intermediate manometers, as 
shown, indicated pressure differences across various longi- 
tudinal sections of the column of sand. 

The gas-liquid mixture discharging from the flow tube was 
conducted through a refrigeration unit to condense it. This unit 
was a thermostatically controlled water bath through which 
the discharge line passed in the form of a coil of copper tub- 
ing. The propane that departed from the cold bath was entirely 
liquid and remained in this state until it returned to the 
column of sand. The total weight rate of flow was determined 
by metering the condensed propane with a rotameter (com- 
mercially known as a Flowrater). This instrument 
brated with liquid propane as the test fluid. over a range of 
Accordingly. the pressure and 


was cali- 


pressures and temperatures. 
temperature at the meter, as well as the meter reading. were 
recorded during each flow experiment. 

The propane that departed from the Flowrater was passed 
through a porous-meter filter and then through a two-cylinder. 
variable-speed, diaphragm-type. reciprocating pump, to in 
crease its pressure. The pistons of the pump displaced lubri- 
cating oil, which, in turn, caused a reciprocating motion of a 
flexible. corrosion-resistant, metal diaphragm. separating the 
propane from the lubricating oil. 

The propane discharged by the pump was conducted through 
the hot bath (thermostatically controlled) to raise its tempera 
ture and was then conducted through another 
filter before re-entering the column of sand. 

Each of the two accumulators shown (Fig. 2) was 


porous-metal 


made 
from two hollow hemispheres of steel, bolted together with a 
flange joint. The two halves of each accumulator were sep- 
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arated at the flange joint by a flexible, fluid-tight diaphragm, 
the upper half of the unit having been filled with nitrogen 
gas in connection with a tank of nitrogen and the lower half 
having been filled with propane in connection with the pro 
pane-flow line. The pressure regulator and two pressure gauges 
at each tank of nitrogen were used to measure the pressure in 
the nitrogen tank and to control and measure the pressure of 
the nitrogen in the accumulator. The accumulator at the right 
(Fig. 2) was used primarily to dampen pulsations in the pro- 
pane-flow line and thereby facilitate reading of the Flowrater, 
and the accumulator on the left was used to increase and 
decrease the quantity of propane in the propane-flow circuit. 


Materials 


rhe propane used in the two-phase flow experiments con- 
tained 99+ 
collecting 


mol per cent propane. The purity was tested by 
vapor-pressure data at a few temperatures within 
the range of interest. The results were compared to the pub- 
lished data of Sage. Schaafsma, and Lacey’ and the agreement 
was good. 
Monterey beach sand was selected as the 
raw material for preparation of the synthetic sand used in the 
propane-flow experiments. Visual inspection of this beach sand 
disclosed that it was composed predominantly of quartz grains, 
with a moderate proportion of feldspars and a very small per- 


A steam-cleaned 


centage of ferromagnesian minerals. It contained a trace of 
calcareous material (apparently fragments of sea shells), but 
no organic matter was evident. 

This raw stock was used as the source from which a supply 
of each of seven size fractions of sand grains was obtained. 
The four fractions desired, nominally 14-28, 
28-35, 35-48, and 48-65 Tyler mesh, were separated from the 
Monterey beach sand directly and in sufficient quantity with 
the appropriate nest of sieves and a sieve shaker. It was neces- 
sary, however. to grind the raw material in a ball mill before 
the three finer fractions, nominally 65-100, 100-200 and 200- 
325 Tyler mesh. could be recovered in sufficient quantity with 


coarsest size 


the sieves and shaker. 

The four coarsest fractions were “de-dusted” by washing 
them thoroughly with tap water. Following this operation, they 
were treated with concentrated hydrochloric acid to remove 
acid-soluble material and washed on a filter with tap water 
until the filtrate would fail to change the color of blue litmus 
paper. Finally, they were dried by placing them in an oven 
at 350°C for 72 hours. 

The fineness of much of the material from the mill caused 
agglomerates to form on the screens during the sieving opera- 
tions with the three finest fractions. These operations, there- 
fore, were followed with a more refined secondary process of 
size fractionation. In this respect, each of the three products 
obtained from the sieving operations was rendered “clean cut” 
by further separation of water, applying Stokes’ law” of set- 
tling velocities. 

The products of the wet-separation process were treated 
with concentrated hydrochloric acid and then washed thor- 
oughly on a filter with tap water until the filtrate would fail 
to change the color of blue litmus paper (as was done with 
the coarse fractions). The acid-free grains were washed with 
isopropyl alcohol to remove water and then with petroleum 
ether to remove alcohol. After the petroleum ether evaporated 
from the grains, they were dried in the same manner as the 
coarse fractions — by placing them in an oven at 350°C for 
72 hours 

rhe finest of 
passed through a 


the seven end-product sand fractions, which 
74 mi- 
(open- 
hydrometer 


200 Tyler mesh sieve (opening 
325 Tyler mesh sieve 
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analysis” and found to contain particles as small as 24 microns 
in Stokes’ diameter. 

The sand used in the propane-flow experiments was made by 
blending the seven foregoing clean-cut fractions. In terms of 
a Tyler sieve analysis, the end-product blend contained, ré 
spectively, 11.0, 8.7, 11.0, 11.0, 12.1, 22.0, and 24.2 per cent 
of 14-28, 28-35, 35-48, 48-65, 65-100, 100-200. and 200-325 
mesh grains. 

\ testing machine (for determining the tensile and com- 
pressive strengths of engineering materials) was used in the 
sand-packing operations, the flow tube being in a vertical 
position in the machine with the packing being done against 
the upstream end. The sand was introduced in 250-gram in 
crements, each having the foregoing end-product grain size 
distribution. After each increment was put in the tube the sand 
was compressed under a compaction pressure of 1,000 psi 
with simultaneous jarring of the tube. In each instance the 
upper surface of the sand was ruffled slightly, before another 
increment was introduced, so that stratification discontinuities 
would be eliminated. The thermocouple junctions were buried 
in the sand, at selected distances along the tube, during the 
packing operation. The average porosity of the resultant col 
umn of sand was calculated to be 27.3 per cent. 


Experimental Procedure 


After the tube was packed and assembled, the homogeneous 
fluid permeability distribution of the column of 
determined with dry air, by application of standard test pro 
cedure. Over a period of several days this test was repeated 
twice and the permeability distribution was found to be un 
changed, indicating that the column was stable. Moreover. the 
test was repeated again with all-liquid propane as the fluid 
medium, between the third and fourth two-phase flow experi 
ments. The results agreed with those obtained in the dry-air 
flow tests. 

Preparatory to two-phase flow experimentation the system 
(Fig. 2) was twice evacuated and purged with the test grade 
of propane gas. With the system filled with this propane gas 
at about 20 psia and room temperature, action was taken to 
introduce liquid propane. A high-pressure container (conven 
tional “gas bottle”) of the test grade of saturated liquid and 
gaseous propane was inverted, connected to the suction line 
on the pump, and elevated so that the saturated liquid propane 
would tend to gravitate into the system. Precautions were 
taken to remove air from the connecting line and a filter was 
assembled in this line to preclude entrance of foreign matter 
into the system. The needle valve (Fig. 2) at the downstream 
end of the flow tube then was closed and the pump used to 
transfer propane from the supply tank to the system. This 
operation was performed at room temperature and as the fluid 


sand was 


pressure increased in the system, gaseous propane condensed 
The maximum pressure was controlled by shutting off the 
pump and by appropriately actuating the accumulator dia 
phragms. This propane-transfer operation was completed when 
the propane-supply tank was disconnected from the pump and 
the system filled at room temperature with liquid propane 
with the accumulator diaphragms about midway between the 
upper and lower extremities of their vertical paths of travel. 

The second operation performed in making the apparatus 
ready for two-phase flow experimentation, was to circulate the 
propane in the system, and simultaneously, to establish a pre 
determined air-temperature distribution in the duct enveloping 
the flow tube, to adjust the temperatures of the hot and cold 
baths as required, to withdraw propane from the propane 
flow circuit, as necessary, with the accumulators, and to regu 
late the pump speed. In this operation, establishment of a 
predetermined pressure and temperature for the liquid pro 
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pane entering the flow tube and establishment of a predeter- 
mined pressure or temperature for the gas-liquid mixture 
leaving the flow tube, fixed the necessary bath temperatures, 
the amount of propane to be withdrawn by the accumulator, 
and the pump speed. 

\s a first approximation a linear temperature distribution 
was established in the air duct, in the direction of propane 
flow. After steady-flow conditions were attained, temperature 
differences, between the air outside the flow tube and the pro- 
pane inside, were measured. The air temperatures then were 
adjusted to match the propane temperatures. When steady- 
flow conditions again were attained, this temperature measure- 
ment and adjustment process was repeated. It was found that 
the inside and outside fluid-temperature distributions would 
remain coincident after two or three such adjustments and 
that little additional regulation of the pump speed and hot 
and cold bath temperatures, was necessary. 

Finally, under conditions of steady flow, with no appreciable 
transfer of heat across the walls of the flow tube and with 
liquid propane entering the tube and a gas-liquid mixture of 
propane discharging from the tube, temperature and pressure 
distributions in the propane, in the direction of flow, were 
measured, as well as the total weight rate of flow of the pro- 
pane. At the conclusion of each of the five experiments the 
propane in the system was vented to the atmosphere, with 
nitrogen from one of the tanks of nitrogen (Fig. 2). 


EXPERIMENTAL AND THEORETICAL 
RESULTS 


Pressure-Temperature Data and the 
Vapor-Pressure Chart 

Fig. 3 depicts a segment of the normal (plane-surface) 
vapor-pressure curve for propane, plotted from data published 
by Sage, Schaafsma, and Lacey.’ Experimentally obtained pres- 
sure and temperature data from the present two-phase propane 
flow experiments are super-imposed on this chart. Each experi- 
mental point shown represents the pressure and temperature 
measured at some transverse section in the sand-packed flow 
tube during two-phase flow experimentation. 

As described previously, the flow would be essentially iso- 
thermal in the all-liquid flow portion of the system. All-liquid 
flow is reflected on the chart for each of the five experiments, 
by the points that are above the vapor-pressure curve and 
form a straight line parallel to the pressure axis. 

It has been indicated further that gas-liquid capillary pres- 
sures and capillary-pressure gradients would be of negligible 
magnitude if experimentally measured liquid pressure-temper- 
ature data, pertaining to the two-phase flow portion of the 
tube, fell on the vapor-pressure curve. As the liquid propane 
was the wetting phase in the experiments, the experimental 
points on Fig. 3 represent liquid pressure-temperature data. 
With the exception of one or two points, data pertaining to 
the two-phase flow portion of the tube fall, within the limits 
of experimental accuracy, on the vapor-pressure curve. Ac- 
cordingly, this method of analysis indicates that the magni- 
tudes of capillary pressures and capillary-pressure gradients, 


are negligible. 


Selection of Relative-Permeability Relations 


Theoretical calculations of pressure distributions and weight 
rates of flow. as made with Equation (22), for comparison to 
experimental results, require that k,/k and ke/k be known 
as functions of pressure. These functions, evidently, may be 
determined indirectly from the proper independent relation- 
ship between k,./k, and k,/k, considering that Equation (24) 
yields k,/k, as a function of pressure. Validity of the theory 
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FIG. 3 — VAPOR-PRESSURE CURVE FOR PROPANE WITH SUPERIMPOSED 
PRESSURE-TEMPERATURE DATA FROM PROPANE-FLOW EXPERIMENTS 


must be assumed, however, if the necessary relationship be- 
tween k, /k,, and k,./k is to be established with the test column 
of sand and two-phase fluid (propane), under conditions of 
adiabatic flow. Some other more appropriate course of action, 
therefore, is required to establish this relationship. 

Based on the findings of Wyckoff and Botset,’ the relative 
permeabilities of unconsolidated sands to either the wetting 
or nonwetting fluid depend almost entirely on the wetting- 
phase content of the sands. In this respect, those investigators 
experimented with a number of two-phase fluids and uncon- 
solidated sands (ranging in homogeneous-fluid permeability 
from 17.8 to 262 darcys), and found that the relative perme- 
abilities were practically independent of the test two-phase 
fluid and were influenced very little by the textural charac- 
teristics of the sands. When &,,/& and k, 
experiments were plotted, on the same chart, against wetting- 
phase content, it was found that the plotted points, as a whole, 


k data from all their 


indicated only one wetting-phase and one nonwetting phase 
relative-permeability curve, the scattering of the points about 
these curves being moderate and probably within the limits 
of experimental accuracy. 
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Consequently, it appears that the Wyckoff-Botset results 
could be used to plot a relationship between k,/k, and k,/k 
that would describe reasonably well the relative-permeability 
characteristics of all unconsolidated sands. One could hardly 
expect this relationship to be exact, however, in describing the 
column of sand used in the propane-flow experiments as it had 
a homogeneous-fluid permeability of the order of one darcy. 
The foregoing findings, nevertheless, disclose that the rela- 
tionship in question should represent a curve that lies in the 
vicinity of the Wyckoff-Botset curve. 

Pursuant to this line of reasoning, data from one of the 
five propane-flow experiments were considered. The inlet tem- 
perature and pressure and the outlet pressure measured in 
this experiment (No. IL). were used with the average curve 
through the Wyckoff-Botset results (Fig. 4) to calculate by 
Equation (22) the theoretical weight rate of flow and pressure 
distribution. The calculated results were compared to the 
experimental results and the curve shifted accordingly. The 
new curve was applied similarly and the curve-shifting proce- 
dure repeated until good agreement of calculated and experi- 
mental results was obtained. The final curve, located in this 
manner, is the propane-flow curve of Fig. 4. It falls well 
within the limits of experimental accuracy of the Wyckoff- 
Botset results. This propane-flow curve then was used in con- 
junction with Equation (22) to calculate total weight rates of 
flow and pressure distributions pertaining to the other four 
experiments. 

In this conection the greatest experimentally-measured total 
weight rate of flow was approximately 100 per cent greater 
than the smallest total weight rate of flow, and the pressure 
and temperature distributions along the test column of sand 
varied widely among the five experiments. Thus, agreement of 
calculated and theoretical results for the four experiments 
used in testing the theory, would indicate that the propane- 
flow curve of Fig. 4 is independent of weight rates of flow 
and fluid properties, and, therefore, dependent only on the 
porous medium. 


Substantiation of Theory 


In making theoretical calculations the pressure and temper- 
ature of the all-liquid stream entering the column of sand 
and the pressure of the gas-liquid mixture discharging from 
the column, were considered as given. The problem was to 
find the corresponding pressure distribution and total weight 
rate of fluid flow. 

In this respect. Table I presents a comparison of experi- 
mentally measured and theoretically calculated pressure dis- 
tributions and total weight rates of flow. Substantiation of the 
theory is based on results pertaining to experiments I, III, IV 
and V, as data collected in experiment II were used as an aid 
in determining the propane-flow curve of Fig. 4. 

Omitting therefore. it is found that for 
experiment I the theoretical total weight rate of flow is 4.8 
per cent greater than the corresponding experimental rate. 
and for experiments II], IV and V. the theoretical rates, 
respectively, are 5.4, 2.5 and 0.7 per cent less than the experi- 


experiment II, 


mental rates. Similarly, the maximum difference between the- 
oretical and experimental pressures for experiments I, IIT, TV 
and V are. respectively. 0.6, 0.6, 2.1 and 0.8 psi. In these 
experiments. in the same order, the total experimentally meas- 
ured pressure differentials between the upstream and down- 
stream ends of the sand column, were 67.9. 37.2, 74.6 and 
58.8 psi. 

Limitations of the experimental equipment and method of 
calculation disclose that the foregoing results (Table I) serve 
as a substantiation of the theory. Regarding the equipment. 


the test pressure gauge at the upstream end of the column of 
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sand was sensitive to one psi, the calibrations of the Flowrater 
indicated that that instrument could introduce errors up to 
about three or four per cent in the total weight rates of flow, 
and calibrations of the thermocouple circuit indicated that 
errors of the order of 1°F were likely to occur in temperature 
measurements. In connection with the method of calculation 
used, experiment Il (Table 1) shows a theoretical weight rate 
of flow 0.3 per cent greater than the experimental rate and a 
maximum difference of 1.5 psi between theoretical and experi 
mental pressures, the total drop in pressure in this experiment 
having been 81.5 psi between the upstream and downstream 
ends of the column of sand. Thus, as data from experiment II 
were used in finding the propane-flow curve of Fig. 4, differ 
ences between theoretical and experimental results for the 
other four experiments could be expected to be as great as 
those cited for experiment IL. Homogeneous-fluid permeabili 
ties measured for various increments of length along the sand 
column, were assumed, in the calculations pertaining to two 
phase flow, to apply to the centers of the increments; hence 
unknown local variations in homogeneous-fluid permeability 
would introduce errors. Finally, graphical methods of analysis 
in the calculations. It is estimated that thess 


two 


were necessary 
methods could introduce errors of the order of one or 
per cent or more. Many of the viscosity data were obtained by 
extrapolation of curves. 

Experimentally determined pressure distributions in the two 
phase flow portion of the column of sand, are represented 
graphically in Fig. 5. This figure shows the ranges of pressur¢ 
used in the experiments, indicates comparative pressure gradi 
ents and delineates the part of the column in which the fluid 
existed in two phases. The lefthand terminus of each of the 
curves defines the transverse section at which vaporization 
method of calculation described 


started, determined by a 


subsequently. 


Details of an Application of Method 
Of Theoretical Calculation 


The method used in making theoretical calculations is pre 
sented by developing theoretical results pertaining to experi 
ment V, which was representative of all the experiments. 

The pressure and temperature of the all-liquid propane 
stream entering the column of sand in this experiment were 
193.2 psia and 96°F, respectively. With these data it is found 
that the enthalpy of the compressed entering liquid is 24.4 
Btu per lb of fluid flowing; thus, the constant enthalpy (7) of 
the process established (see Fig. 1). 
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FIG. 4— RATIO OF RELATIVE PERMEABILITIES AS FUNCTIONS OF REL 
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Table I—Comparison of Theoretical and Experimental 


Pressure Distributions and Total Weight Rates of Flow 


Experiments 
il il IV 
Experimental and theoretical pressures, psia 


Exp. Theor. Exp. Theor. Expt. Theor. Exp. Theor. Exp. Theor. 
193.2 
189.0 
184.9 
179.7 
169.4 
164.2 
145.9 
140.3 
134.4 


215.6 
210.4 
205.4 
199.1 
186.9 


193.2 
188.5 
184.8 
179.6 
170.2 
154.7 
146.4 
140.4 
134.4 


0.000 210.8 210.8 22 2 190.3 215.6 


206.0 205.8 2 2.3 7.3 187.6 210.3 


200.9 216 . 85. 205. 
194.9 209.7 209.9 7 b «6198.8 
196.6 73 3.2 186.2 
175.7 5 a4. 168. 167.5 
57 163.8 59. 59.2 158. 156.7 
150.1 155.3 562 149. 149.1 
142.9 146.3 146.3 53. 53. 141.0 141.0 


Experimenta! and theoretical total weight rates of flow, lb/hr 
12.45 1435 14.40 7.045 6.664 13.42 13.08 10.36 10.29 


11.88 

For saturated liquid propane an enthalpy of 24.4 Btu per 

lb corresponds to a pressure of 180.1 psia. At this pressure, 
therefore, vaporization started in the column (see Fig. 3). 


To find the transverse section in the column, at which the 
pressure was 180.] psia, it is noted first (Fig. 3) that the 
flow was essentially isothermal in the part of the column up- 
stream from the place of initial vaporization. Moreover, the 
measured total weight rate of flow was 10.36 lb per hr, and 
at distances of 1.669 and 2.510 ft along the column, the meas- 
ured pressures were 184.8 and 179.6 psia, respectively. In the 
indicated incremental length of sand column, the bulk cross- 
sectional area was 0.0237 sq ft and the homogeneous-fluid 
permeability 1.001 x 10" sq ft. Graphing of published data*’ 
discloses that the product of specific volume and absolute 
viscosity was, 62.18x 10° ft+sec for the liquid propane that 
was flowing between the transverse sections at which the pres- 
sures were 184.8 and 180.1 psia, respectively. Hence, applica- 
tion of Darey’s law in the conventional manner to homogene- 
ous-liquid flow, gives 2.440 ft from the upstream end of the 
sand column as the distance to the transverse section where 
vaporization started (see Fig. 5). 


Figs. 6 and 7 are, respectively, the pressure-function and 
integrated-pressure-function charts. These charts are necessary 
in determining the theoretical total weight rate of flow and 
pressure distribution. Data used in the preparation of these 
figures are entered in Table IL and were developed as dis- 


( ussed below. 


In Table Ll, the pressures at the upstream and downstream 
ends of the column, 193.2 and 134.4 psia, respectively, were 
measured; and, the pressure at which vaporization started 
(180.1 psia) was determined by the foregoing procedure. The 
pressure range from 180.1 to 134.4 psia was subdivided arbi- 
trarily as shown in Column 1 so that data could be developed 
for constructing Figs. 6 and 7. 

In the all-liquid flow portion of the sand column k,/k ob- 
viously would have been zero and k,,/k unity; hence, for this 
part of the system F(P) would be the reciprocal of 4,7, (see 
Equations 2] and 22). This reciprocal was calculated to be 
16.08 x 10° ft’ . and the horizontal segment of the F/(P) 
vs P curve (Fig. 6) thereby became established. 


* sec 


lable Il contains two entries of F(P) for the pressure 180.1 
psia. At the upstream side of the transverse section at which 
vaporization started, the fluid in the sand was entirely in the 
liquid state, and therefore, F(/P) for that section would be 
16.08 x 10° ft''+ see’ if it were approached from the upstream 
side. In accordance with the concept of an “equilibrium gas 
saturation,” however, the pore space in the sand at the section 
in question would have contained a fixed quantity of static 
gas but no moving gas. This means that k,/k would have been 
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FIG. 5 — EXPERIMENTALLY-DETERMINED PRESSURE DISTRIBUTIONS IN 
TWO-PHASE FLOW PORTION OF TEST COLUMN OF SAND. 


zero and k,/k less than unity at the section. Hence, if the rela- 
tive permeability of the liquid at the section is denoted by 
(k,/k).. “the equilibrium permeability,” and the section ap- 
proached from the downstream side F(P) for the section would 
be (hk, /k)./(m.%,).* 

For unconsolidated sands, Wyckoff and Botset found that 
(k,/k), t= 0.79. This value was used in the present analysis. 
\ considerable error in this value would have no appreciable 
effect on the theoretically calculated results (Fig. 6). 

With respect to the two-phase flow in the column of sand, 
published data*’ on saturated liquid and gaseous propane 
were used to plot (#:¥t/meva) and (H-H,)/(H,.-—H,) as 
functions of pressure. These plots were the source for data 
entered in Columns 2 and 3 of Table II. This information and 
application of Equation (24) led to the calculated values of 
k./ke shown in Column 3. In turn, the data in Column 3 were 
used with the propane-flow curve of Fig. 4 to determine the 
corresponding relative permeabilities (k,/k and ke/k). Fur- 
ther use of the foregoing published data to plot 4, ¥,, and “ei. 
as functions of pressure, together with knowledge of k,/k 
and k,/k as functions of pressure, resulted in the entries in 
Columns 5 and 6. Corresponding entries in these columns were 
added to give the F(P) entries in Column 7, which were used 
to construct Fig. 6. 

Data for the construction of the integrated-pressure-function 
chart (Fig. 7) are entered in Column 7 of Table Il. These 
data were obtained by dividing the total area under the curve 
in Fig. 6 into vertical strips and determining the area of each 
strip by application of either Simpson's rule or the trapezoidal 
rule. As the widths of the strips were governed by the pressure 
entries in Column 1, knowledge of their areas led directly to 


*It is of interest to find that there would be a discontinuity in the fluid 
pressure gradient at the transverse section under consideration. Applying 
Darcy's law, the pressure gradient at the upstream side would be 

(Wyuzy¥,/Ak) and at the downstream side ~(Wy,7,)/[Ak(k,/k) .1; 
hence the ratio of upstream to downstream gradients would be (k, /k) 
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the entries in Column 7, which were used to construct Fig. 7. 


To calculate the theoretical total weight rate of flow, 


(22) was put in the form 


P, 
fh 

if F(P)dP. . . (29) 
Ax 

P, 


in which subscript i refers to the eight length increments into 
which the column was divided;* P, and Py are the fluid pres- 
sures respectively at the upstream and downstream ends of 


Equation 


e 
o 


i ] 


the column. 


The integral term in Equation (29) represents the total area 
under the curve in Fig. 6 and was calculated to be 706.42 x 10° 
Ib + ft (see Table II, Column 8). The term in brackets 
was calculated to be 2.473x 10” ft*. Hence, the theoretical 
total weight rate of flow was found to be 2.857x 10° Ib per 
sec or 10.29 lb per hr as compared to 10.36 lb per hr for the 


* sec 


measured rate. 


To find the theoretical pressure at the upstream end of the 
n-th increment (n = 1, 2, 3, 4, 5, 6, 7 or 8), the upper limits 
of the integral and summation in Equation (29) are changed 
to P,. Under these conditions, the value of the integral may 
be calculated, as the theoretical weight rate of flow will be 
known. The calculated value of the integral, then, is used in 
conjunction with Fig. 7 to find P,, the theoretical pressure at 
the upstream end of the n-th Table I 
results determined in this manner. 


increment. contains 


*These eight incremental lengths were limited, respectively, to the fol- 
lowing ranges of distance (/\x) along the sand column (measured in ft 
from the inlet face): 0.000 to 0.833, 0.833 to 1.669, 1.669 to 2.510, 2.510 
to 3.351, 3.351 to 4.187, 4.187 to 4.570, 4.570 to 4.778, and 4.778 to 5.028. 
Starting at the upstream end of the column and proceeding downstream, 
the homogeneous-fluid permeabilities (k) of these increments were, 
respectively, 1.191, 1.244, 1.001, 0.8950, 0.8812, 0.9110, 0.7951 and 1.006 
sq ftx 10 Considered in the same order their average bulk cross-sec- 
tional areas (A) were, respectively, 0.02042, 0.02036, 0.02037, 0.02044, 


0.02025, 0.02027, 0.02027, and 0.02027 sq ft 


o @ re) rN) rs o 


PRESSURE FUNCTION, F(P), ff''sec'x 10 © 
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PRESSURE, P, psia 
FIG. 6 — GRAPHIC REPRESENTATION OF PRESSURE FUNCTION (BASED 
ON THEORETICAL CALCULATIONS PERTAINING TO EXPERIMENT V). 
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CONCLUSION 


this paper 
namic aspects of the steady flow of a single-component gas 


liquid mixture, through porous media. The underlying. funda 


The subject matter of concerns the thermody 


mental two-phase fluid-flow theory pertaining to this subject i 
developed, and laboratory experiments substantiating the the 
ory are described. 

It is pointed out in the paper that the thermodynamic pro 
ess followed by fluids while 
porous underground rock formations is of significance in pre 


vil-reservoir moving through 
dictions of oil-reservoir behavior. Moreover, the flow of multi 
hydrocarbon fluids in these rocks is indicated to 
be intermediate isothermal and adiabatic, but not 
susceptible to exact classification owing to a lack of necessary 
detailed, heat-transfer data on rocks fluid 


component 
between 
subsurface and 
mixtures. 

In an effort to perform the described experiments with the 
simplest appropriate two-phase flow system, it was decided that 
a single-component fluid should be used. This decision mad 
isothermal flow impossible; hence. on the basis of the fore 
going considerations, adiabatic flow was adopted for study 

In this respect, the experiments reported relate specifically 
to the steady adiabatic flow of propane in two phases throug! 
a horizontal column of sand, mounted in a specially designed 
apparatus. The propane, in all the experiments, entered the 
column entirely as a liquid, started to vaporize an appreciable 
distance downstream from the inlet, and then flowed in two 
phases throughout the remainder of the column. 

Regarding the theory, the problem posed was to calculat 
the total weight rate of flow and the pressure distribution along 
the column of sand, for a given pressure and temperature of 
the all-liquid entering stream, a given pressure or temperature 
of the discharging gas-liquid mixture. and given sand and 
fluid characteristics. A comparison of theoretical and experi 
mental results has been presented and discussed to show that 
the theory may be applied to solve this problem. 

The theory, as developed herein. is premised on work done 
by O’Brien and Putnam.’ The present paper shows, however 
that equations developed by those investigators for calculating 
wetting-phase content and used by them in determining weight 
rates of flow and pressure distributions were not correct 


Table II 
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FIS. 7 — GRAPHIC REPRESENTATION OF INTEGRATED PRESSURE FUNC- 
TION (BASED ON THEORETICAL CALCULATIONS PERTAINING TO EX- 
PERIMENT V 


\ttention is given in the present paper, moveover, to gas- 


I:quid interfacial phenomena. The analyses made in this con- 
iection indicate that capillary pressures and their gradients 
are ot negligible magnitude for single-component two-phase 
flow systems of the kind investigated. 


Furthermore, the agreement among theoretical and experi- 
mental results indicates that the assumption of phase equilib- 
rium, made in the development of the theory, is valid for the 
system studied 

In conclusion, the material offered in this paper should lead 
to a better understanding of oil-reservoir behavior in that the 
nvestigation pertains to the flow of a two-phase hydrocarbon 
flrid in sand. Substantiation of the theory suggests that it 
would be desirable to follow the reported work with a similar 
tudy of the two-phase flow of a multi-component hydrocarbon 
gas-liquid mixture through sand. 


Data for Construction of Pressure-Function and Integrated-Pressure-Function Charts (Figs. 6 and 7, 


Respect iv ely ) 


(4 


} 
J 


dimension- dimension- dimension 


less less less 


0.3622 
0.3631 
0.3704 
0.3779 
0.3850 
0.3916 
0.3978 
0.4037 
0.409] 
0.4137 


0.1067 3.032 
0.1047 3.105 
0.09386 576 
0.08224 217 
0.07038 086 
0.05908 238 
0.04762 7.956 
0.03598 82 
0.02417 : 
0.01293 

180.1 0.00000 
180.1 7 
193.2 


*(k,/k)o = 
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0.79 — from Wyckoff - Botset 


results. 


(8) 
pP 
F(P) { 
P, 
se Ib 
x 10° 


F(P)dP 


ft sec ' 
x 10° 


0.4154 
0.4174 
0.4011 
0.3784 
0.3489 
0.3181 


706 
0.2788 


3.892 
3.987 
4.273 
1.601 
1.959 
5.385 
5.854 
6.481 
7.329 
8.679 
12.71 


16.08 


0.00 
3.57 
33.31 
65.26 
99.68 
136.92 
177.38 
221.79 
271.34 
328.57 
103.08 
103.08 
706.42 


0.233 

0.1771 

0.1122 
12.71% 0.0000 
16.08 ’ 
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DISCUSSION 


By E. C. Babson, Union Oil Co. of California, Member 

Miller has 
between experimental results and analytical theory. indicating 
that rapidly under the 
conditions of his experiments. 


11ME 


succeeded in obtaining remarkable agreement 
phase equilibrium is attained very 


rhe relative permeability data used in calculating the pres- 


sure gradients, however, were not established for the sand 


column used in the experiments but were obtained by shifting 
K./K 


rhe relation established in this manner was used in the 


Botset’s curve to reach an overall fit on one of the 
runs, 
calculations for all runs. 

rhere is, therefore. a possibility that the agreement between 
theory and experiment is due to compensating errors. When 
consideration is given to the excellent agreement at each pres- 
sure tap for each run. it must be admitted that this possibility 


is remote, but it still exists and should be recognized. 


AUTHOR'S REPLY TO MR. BABSON 


Babson’s comment regarding relative-permeability data is 
well taken 
more detail in the paper. 
includes a statement to the effect that 


Perhaps this matter should have been discussed in 
rhe paper validity 
of the theory must be assumed if the relationship between 
and kik is to be established with (and for) the test 
column of sand. under conditions of adiabatic flow with two- 
The Wyekoff-Botset unconsolidated-sand data 
zed in the described manner, therefore, as a recourse. 


k k 
pha e propane 
were ut'l 

If this procedure had not been followed, then the collection 
of independent relative-permeability data for the test sand 
column would have been necessary. Considering the high sen- 
sitivity of calculated weight rates of flow and pressure distribu- 


tions to such data. it is evident that an independently-deter- 


mined kk, (ky vs Aish would 


-hifting (within the limits of experimental accuracy) before 


curve most probably require 
an acceptable agreement of theory and experiment were ob- 
tained. As this shifted curve would have about the same the- 
eretical significance as the propane-flow curve of Fig. 4, it too, 
would give rise to a remote possibility of an unsound agree- 
Accordingly, one concludes 
from the collection of 

oe 2: & 


ment of theory and experiment. 


that no important jain would result 


independent relative-permeability data. 
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CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA, CALIF. 


INTRODUCTION 


The transient flow of gases in long pipe lines is a problem 
of industrial interest. The present discussion deals with the 
application of the conservation of momentum and material to 
the transient flow of gas in conduits which are long with re 
spect to their diameter. Solutions of the resulting equations 
include consideration of variations in friction and 
volume of the fluid as functions of time and position along 
the conduit. The changes in state along the pipe line are 
established as a function of time. The methods proposed may 
be extended to other situations. 


spec ifi 


The flow of compressible fluids under steady conditions 
has been considered in detail. Primary emphasis was placed 
on the development of simple means of taking into account the 
non-uniform flow’ realized in long conduits under such condi 
tions. Binder’ discussed the essential differences between the 
flow of compressible and noncompressible fluids. Bonilla 
proposed direct methods of solution for equations describing 
isothermal flow in long pipes. Similar considerations limited 
to steady flow have been discussed by Ruth.’ Palsgrove’ sug 
gested a graphical means of solving specific problems relating 
to the flow of gases in conduits of uniform section. Some of 
the basic concepts of the flow of a compressible viscous fluid 
in straight tubes have been described by Kuo.’ Problems of 
this type require knowledge of the equation of state of the 
fluid. In 


the fluid follows the relationship ascribed to a perfect gas 


many cases it has been convenient to assume that 

The treatment of unsteady flow of compressible fluids is less 
complete. This state of affairs results from the fact that the 
differential equations relating pressure and flow rate with time 
and position for unsteady flew are nonlinear in character 
Muskat and Wyckoff presented an excellent treatment of the 
flow where inertia 
Somewhat 


forces were not of 
controlling importance. earlier Hurst” 
the unsteady flow of fluids in petroleum reservoirs. Hethering 
ton, MacRoberts, and Huntington” established experimentally 
the pressure distribution and accumulative flow rate of natural 


situation in laminar 


dis« ussed 


gas through an unconsolidated sand. Recently Joffe con 
sidered the specific problem of the flow of natural gas in rela 
tively long pipe lines for unsteady conditions with emphasis 
upon the economic feasibility of the storage of gas in such a 
reservoir. These discussions were based upon approximations 
of the partial differential equations applying to the conditions 
of flow. The results were descriptive of the physical situation 


evelic 


and indicated the effectiveness of large pipe lines a 


References given at end of paper 
Manuscript received by the Petroleum Branch Oct. 27, 1950 
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storage facilities for natural gas. The advent of analog com- 
has increased materially the practicality of solving 
partial differential Several iterative 
numerical and graphical methods’*’’” have been developed to 
solve nonlinear equations. However, it appears that the pair 
of simultaneous partial differential equations in- 
volved in problems of this nature are not directly susceptible 
to solution by presently available analog computer or formal 
iterative numerical techniques. Nevertheless, graphical analy- 
ses afford a means of obtaining data of engineering value 
without special equipment. 


puters’ 


nonlinear expressions, 


nonlinear 


MATHEMATICAL CONSIDERATIONS 


Since the length of the conduit is large in comparison to 
the cross section in the cases to be considered, it is possible 
to treat the characteristics of the flow as independent of posi- 
tion in the cross section of the conduit. This possibility follows 
since the energy connected with the velocity profile is small 
compared to that required to overcome the friction associated 
with its movement throughout the length of the line. Average 
kinetic energy, momentum, 
and specific volume may be assigned to the fluid at any par- 


values of the pressure, velocity, 
ticular cross section. From a mathematical viewpoint this may 
be expressed in the following way for velocity: 

SS BPRS 5 6: Soa ee eo ee 
The uncertainties introduced by such assumptions usually are 
negligible when the ratio of length to diameter is greater than 
1,000. On this 
be considered to occupy the entire cross section of the conduit. 


the order of basis a volume element may 


Such an element is presented in Fig. 1. 


























FIG. 1 — VOLUME ELEMENT IN A LONG PIPE LINE 
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\ material balance applied to the volume element yields 


the following for a conduit of cross-sectional 


area A: 


expression 


ce l ojn 
Al 1 AX 


The foregoing equation represents the fact that the differ- 


0 er ee eae oe ks (2) 


ence in the material entering and leaving a particular part of 
the pipe line in a given time is equal to the change in the 
quantity of material within the section of the pipe line under 
consideration. 

A force balance upon the differential element of volume 
shown in Fig. 1, which is equivalent to assuming conservation 
of momentum, is given by the following equation: 


l am eC m rd 
' (v pe ae 
zi 00 ax goA { 


Following conventional practice.” the shear per unit length is 
empirically related to the conditions of flow and the friction 
factor in the following manner: 


» 
{o) 


rrd i m 
— ee 


1 2og dA 


\ combination of Equations (3) and (4) vields the expression 


‘ m fan 
P 0 
gi 36 ox goA 2egdA 


For the majority of conditions encountered in the transient 


l em 


flow of natural gas in long pipe lines the imertia effects are 
small. and it is believed that except in special circumstances 
where numerous large changes in velocity are encountered the 
forces associated with the acceleration of the flow can be neg 
lected without introducing appreciable error. In the examples 
to be considered the neglect of these inertia terms did not 
introduce errors greater than 0.8 per cent in the reported 
values of pressure. If the inertia effects are neglected, Equa 
tien (5) simplifies to the following form: 


eP 
ax 2ogdA 


The evaluation of Equations (2) and (6) requires a knowl 
1 1 


fin 
ve 2 : . (6) 


edge of the relationship between pressure and specific weight 
of the gas. In the present instance it is not necessary to em 
ploy a complex equation of state and the following expression 
ippears adequate for the present discussion: 

P 


bT— PB 





22 MILES 








25.13 INCH DIAMETER 


EXAMPLE A 





165 MILES 








24.88 INCH DIAMETER 
, 


+ 24.88 INCH DIAMETER 


EXAMPLE 8B 











FIG. 2—SCHEMATIC DIAGRAM OF PHYSICAL SITUATION FOR EXAMPLES 
A AND B 
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FIG. 3— PRESSURE AS A FUNCTION OF POSITION UNDER STEADY 
STATE CONDITIONS IN EXAMPLE A 








rc 


Or. neglecting terms of order this may also be 


(PB) bT) 
form 


l B 
(r+ fr) 
hT bT 


dual volume V has been considered a constant for the 


written in the 


(8) 


The re 
range of of interest as is indicated in the third term 


ind Equatien (8). 


pressure 
of Equation (7 


Description of the Situation 


Example A’ Example B 
Length of pipe line 
mile- 122 

0.6442x 10 
2.094 

60 

815 

300x10°** 


214* 
1.131x10 
2.072 


60 


leet 
Diameter of conduit. feet 
len perature degrees I 
l pstream pressure 
Flow rate. cu ft/24 hr 


( omposition of gas 


ps! 
300x 10 
mole fraction 

0.94 

0.04 

0.02 
17.28 
0.2799x 10 
0.05471 


methane 
ethane 
propane and heavier 
19.79 
0.2108x10 
0.0495 


ular weight 
Viscosity of gas. Ib. sec 
Residual volume. cu ft/lb 

Constants of Equations {9) and (10) 
0.019724 
0.07534 
0.1205 


Average mole 


sq ft 
specific 


0.11754 
0.1454 


with the 


cross section a 


length 
ame 


of the 


final 49 miles made up of a pair 
the single conduit making 
vstem 
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SPECIAL CONDITIONS 


The remaining portion of this discussion is concerned with 
two special problems which differ materially as to the physical 
situations. The units of measure used in the two examples 
have been chosen somewhat differently. However, the force 
length-time system of dimensions has been employed. The sit 
uation associated with each of the pipe lines is presented in a 
part of Table I, and in Fig. 2 is depicted a diagrammatic 
sketch of their arrangement. 

The friction associated with the flow was established from 
the correlation of Pigott.” which is based in part upon the 
experimental work of Nikuradse.” The friction factor 
pirically related to the Re ynolds number in the following way 

f=eceRe* . ~ « (9) 
It is assumed that viscosity is inde >pe dent of pressure over 
the range of conditions encountered. If the diameter of the 
line is uniform, it follows that the friction factor becomes a 
single-valued function of the weight rate of flow 


cated below. 


is em 


as is indi 


‘=e; - ee 
Further consideration of the proble m necessitates treatment of 
a specific set of boundary conditions and the discussion will 
be applied to each situation separately. 


Example A 

The physical situation involved in this problem is set forth 
in a part of Table I. It is desired to determine the variation 
of pressure with respect to time at the exit of a pipe line 122 
miles in length when flow suddenly is curtailed at the inlet 
while gas is withdrawn continuously at the outlet at a constant 
rate of 300 x 10° cu ft per day.* 

The viscosity of the gases employed is reported in a part of 


All volumetric measurements relating to the quantity of gas correspond 
to a temperature of 60°F and a pressure of 14.73 psia. The indicated 
volumetric rate corresponds to 160.7 Ib of gas per second. 
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PRESSURE DISTRIBUTION FOR UNSTEADY CONDITIONS IN 
EXAMPLE A 


FIG. 4 
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INCH 


LENGTH FEET x 10% = 








PRESSURE POUNDS PER SQUARE 

















100 200 
TIME MINUTES 


FIG. 5 — VARIATION IN PRESSURE WITH TIME FOR EXAMPLE A. 


"able | and is based upon experimental measurements of vis- 
In Example A the friction factor is 
related to the weight rate of flow by the following equation: 

f = 0.019724m°"” (11) 
The associated values of the constants of Eeustions (9) and 
(10) are shown in a part of Table I. In this instance the spe- 
cific volume at 60°F may be expressed in terms of the pressure 
by the 


cosity of similar gases.” 


relation: 
l 322.72 
J J ) i re | 
pP 


o 

and (10) followed by 

appropriate integration results in the following expression for 
steady flow: 

P,-P b7 bT/B-P. am?" 

B BOT /B-P, 

In the a pertect 

the assumed constant B 


\ combination of Equations (6). (7). 


(Z,~2Z,) . (13) 


2edA* 
gas where the residual volume, i.e., 
Equation (13) reduces to a 


case of 
Is zero, 
simpler form. 


abTm®™ ( X,- a) 
gd A’ 


specific volumetric data of Equation 
(12) in Equation (13) and the substitution of the boundary 
conditions listed in Table I yield the following expression for 
along the pipe line under steady conditions of 
flow as a function of the prevailing pressure: 


1.2044 x 10 f 
15 — 
(0.019724) in ' 5898.7 


P? -P. (14) 


Ibe application of the 


the distance 


y¥-X 
-PA 
: (15) 
In Fie. pressure as a func tion of position along 
the pipe line for a flow rate of 300 x 10° cu ft per 24 hr. 

In the case of unsteady flow both Equations (2) and (6) 
satished. Under circumstances the pressure 
gradient may be expressed in a slightly different form from 
that orig 


> is shown the 


must be these 


inally employed in Equation (6). 
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VARIATION IN OUTPUT RATE FOR EXAMPLE B 


P h7 ( PB \) a : 
a. P bT 2d 


The change in pressure with respect to time may be obtained 


20 24 


(16) 
directly from Equations (2) and (7). 
aP bT ( PB ') em 
ee { b7 ex 


Both Equations (16) and (17) 
the pipe line. An iterative graphical procedure was employed 


(17) 
are applicable at all points in 


to construct a solation for these equations which satisfied the 
required boundary conditions. The details of a similar graphi 
cal operation are set forth in the latter part of this paper 
The situation at various times after the flow is interrupted 
is shown in Fig. 4 upon the assumption of perfect gas beha 
vior. It is apparent that during the initial period after shutoff 
Fig. 4 


also depicts the pressure distribution after an elapse of 280 


there has been no change in the downstream pressure 


min. At this time the exit pressure bas decreased to 115 psi 
to be the 
for the operation of the transmission line. Fig. 5 shows the 
the pressure 
the pipe 
pressure at the exit is obtained for a period of about 45 min. 


which was considered minimum pressure feasible 


estimated variation in with respect to time at 


various points along line. Litthe or no change in 
Subsequent to that interval, the exit pressure decreases at an 
increasing rate until the effectiveness of the line as a storage 


facility is substantially exhausted after about 280 min 


Example B 
The this 
markedly different from that of the preceding one. The gas 


physical situation considered in example is 
entered a pipe line 214 miles in length which was made up 
of 165 miles of single conduit and 49 miles of two parallel 
conduits of the same cross section as that of the first part of 
the line. This design was employed to increase the storage 
capacity of the system under transient conditions. The system 
was operated at an input rate of 300 x 10° cu ft per 24 hr. The 
output rate of flow varied as indicated in Fig. 6. In the inter 
est of ease of computation the continuous distribution indi 
eated by the dotted curve was simplied to the discontinuous 
the full line. It 


the pressure distribution in the pipe line during a 24-hr period 


relation shown by was desired to determine 
corresponding to the withdrawal rates indicated in Fig. 6 
while the input flow rate remained constant at 300 x 10° cu ft 


per day. The total withdrawal per 24 hr was 300 x 10° eu ft. 
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In this instance the specific weight of the gas at 60°F is 
related to the pressure in the following way: 
0.003547 (P + 0.000176 P*) 


factor 


(18) 
rhe was related to the Reynolds number in 
accordance with Equation (9). 

f = 0.11754 Re*™ ete prea (19) 
\ combination of the foregoing isothermal equation of state 


friction 


and the expression for friction with Equations (2) and (6) 
leads to the following expressions in which the distance is 
the 
feet per day:* 


miles, time in hours, and the flow rate in 


millions of cubic 


expressed in 


ab 
(P + 0.000176 P*) + 33.014 


( 


(P? + 0.0001173 P*) + 0.10788 V"™" = 0 


Z 


Under steady-state conditions Equation (20) 


vanishes and 
Equation (21) may be integrated to yield: 


[P* + 0.0001173 P*]' 0.10788 V(X YX.) 
In Fig 


line for two limiting conditions of steady flow rates. The first 


(22) 
7 is presented the pressure distribution in the pipe 


is associated with the maximum entrance pressure of 970 psi 
which yields an exit pressure of 518 psi. The second condition 
is related to the minimum acceptable exit pressure of 265 psi 
which necessitates an entrance pressure of 870 psi. These con- 
ditions correspond to the limiting range of steady operation 
which permits the transmission of 300x 10° cu ft of gas per 
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FIG. 8 — PRESSURE VARIATION FOR TRANSIENT CONDITIONS IN EX 

AMPLE B. 


24 hr. It was assumed in the last 49 miles of the conduit that 
the flow was divided equally between the two conduits. 

The discontinuity in the slope of the pressure distribution 
curves in Fig. 7 corresponds to the point where the flow di 
vides between the two identical conduits. The transition from 
steady to cyclic exit conditions involved secondary transients 
which were also completely dissipated by the end of the first 
‘ yele. 

The distribution of pressure in the pipe line is depicted in 
Fig. 8 as a function of time and position. At the entrance 
to the conduit a cyclic variation in pressure from 890 to 970 
psi is obtained during a given 24-hr period. At the exit the 
pressure changes from a minimum value of 270 to a maximum 
of 640 psi. The discontinuous nature of the slope of this curve 
results from the assumption shown in Fig. 6 of discontinuous 
changes in the demand rate of flow. In actual practice the rate 
of flow would not change discontinuously and the rate o 
change of pressure with time at the exit would be continuous 
The pressure distribution in the pipe line for several different 
periods in the cycle constitutes Fig. 9. 


A GRAPHICAL METHOD OF SOLUTION 


The pattern of graphical operations used in constructir 
solutions for the examples described in this paper was devised 
only after effort had been expended in attempting to apply 
other approaches. An explanation of the step-wise solution 
employed in connection with Equations (20) and (21) follows 
1. Consider the pipe line at some known conditions through 

out its length. In the present ca-e initial conditions corre 

sponding to steady-state flow of 300 x 10° cu ft per day wer 
chosen. 

Permit the flow to accordance with thie 


output vary in 


prescribed cycle and graphically construct a family of 
curves representing V’ vs X for various arbitrarily chosen 


times. This step appears somewhat difficult. but in practice 
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it may be accomplished crudely by using “intuitive guesses” 
to approximate the behavior of the system. 

3. With the aid of Equation (21) compute P* + 0.0001173 P 
for several values of X. This may be accomplished by 


graphical integration of y’ with respect to \. 

|. Plot curves representing the quantity P + 0.000176 P’ for 
the range in values of Y. Graphical integration of these 
curves with respect to \ will determine the total gas con- 
tent of the pipe line at each of the arbitrarily chosen times 
adopted in Step 2. 

5. Compute from the input and output flow rates the time 
corresponding to each curve in Steps 2 and 4 by correlating 
the cumulative difference in these flow rates with the net 
change in gas content of the pipe line. 

6. Cross plot P + 0.000176 P* with respect to time and differ- 
entiate the resulting curves with respect to time. 

eV’ 


Compute with the aid of Equation (20) and integrate 


graphically plots of this quantity with respect to X so as to 


obtain a new set of relations between flow rate, position, 
and time. 

the process beginning with Step 3 may now be repeated to 
a second approximation. In this example the conver- 
sufficiently rapid so that only three 


iterations of the above-described steps were required to yield 


obtain 


gence was found to be 


a substantial convergence. In this manner, starting from the 


assumed initial conditions, computations for the second ex- 


ample in this paper were carried through two complete cycles, 
corresponding to 18 hr. It was found that the 
pressure distribution at the end of the first cycle did not 


a period ot 


differ significantly from that at the end of the second cycle. 
therefore, that the assumption of an arbi- 
trary the pipe line did not materially 
nfluence the behavior of the system after the first cycle. 


It was concluded 


initial condition for 
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SUMMARY 


The foregoing Examples A and B demonstrate the fashion 
in which the partial differential equations describing the con- 
servation of momentum and material may be applied to the 
prediction of the pressure distribution in long pipe lines carry- 
ing gas under transient conditions. Example A considers a 
problem wherein the flow at the entrance to the pipe line is 
interrupted and the second predicts the cyclic pressure dis- 
tribution to be anticipated in a gas transmission line with a 


constant input and a cyclic demand rate. 
The graphical procedures outlined here for solving these 
problems may be expected to yield satisfactory results for 


wide variations in flow conditions, provided situations involving 
shock waves and supersonic flow do not arise. 
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NOMENCLATURE 


ross-sectional area, sq ft 

onstant defined in Equation (10) 
onstant defined in Equation (7) 
pecific gas constant 

onstant defined in Equation (9) 
liameter of pipe line, ft 

friction factor 

acceleration resulting from gravity. ft per sec per 
sec 

total weight rate of flow, lb per sec 
exponent of Equation (9) 

pressure, psi 

thermodynamic temperature, “R 
velocity, ft per sec 

residual specific volume, cu ft/lb 


volumetric flow rate. cu ft x 10° per 24 hr 
pecific volume, cu ft/lb 
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distance along pipe line. ft 
di-tance along pipe line, miles 
shear at boundary of flowing stream, lb per sq ft 
time, sec 
time, hr 
specific weight, lb per cu ft 
/n —natural logarithm 
log —- Briggs’ logarithm 
Re Reynolds number 
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NOMINEES FOR AIME OFFICES, 1952 


Clayton G. Ball, chairman of the Nominating Committee 
for 1952 Officers, has announced nominations for president, 
president-elect, vice-president and directors for 1952. As 
provided in Article [IX of the by-laws, 25 members or asso- 
ciate members may sign and transmit to the secretary's 
office prior to Sept. 1 “any complete or partial ticket of 
nominees,” should they wish othcr candidates to be con- 
sidered. In such instance, a letter ballot will be forwarded 
to all members in good standing in the United States, Can- 
ada and Mexico, tabulating both the official ballot and any 
supplementary nominations. If no nominations are received, 
no letter ballot will be printed, and nominees on the official 
ballot shall be declared duly elected at the meeting of the 
Board of Directors of the Executive Committee in November. 


FOR PRESIDENT 


MICHAEL L. HAIDER attended the University of South 
ern California and Stanford University, from which he was 
graduated in 1927 with a BA degree in chemistry. He worked 
ior the Pacific Portland Cement Co., Pan-American Petroleum 
Co., Richfield Oil Co., and the Carter Oil Co., before joining 
the Standard Oil Development Co. in 1938 as manager of 
production engineering and research. In 1945 he became 
executive assistant in the producing department of Standard 
Oil Co. (New Jersey). In 1946 he was transferred to Imperia! 
Oil Ltd., Canadian affiliate, with headquarters in Toronto. He 
is currently a vice-president and member of the board of 
directors. Haider was chairman of the Publications Committee 
of the Petroleum Division for several years. He was chairman 
of the Division in 1945, and has held numerous AIME com 
mittee posts. He is now a director of the Institute. 


FOR PRESIDENT-ELECT 
ANDREW FLETCHER, president of St. Joseph Lead Co., 


was born in 1895 in New York, N. Y., to which his pater 
nal grandparents had come from Scotland in 1825. He attended 
Hill School and the Sheffield Scientific School, Yale Univer- 
sity, where he received the PhB degree in mechanical engineer 
ing. In 1949 the University of Missouri conferred upon him 
an honorary degree of Doctor of Engineering. Fletcher has 
been treasurer of AIME since 1944. He is also president of 
Mine La Motte Corp., and Bonne Terre Farming and Cattle 
Co.; vice-president of Missouri-Illinois Railroad Co.; 
director of the American-Hawaiian Steamship Co. 


and 


FOR VICE-PRESIDENTS 


OLIVER B. J. FRASER, assistant manager of the Devel- 
opment and Research Division of the International Nickel 
Co., New York, N. Y., is a past chairman of the New York 
Section and has served on the following AIME committees: 
Ways and Means, Mining and Metallurgy (as chairman) 
President’s Advisory, Nominations for 195] AIME Officers (as 
chairman); Auxiliary Publications Committee of the Metals 
Division (as chairman), and Program Committee, IMD. He is 
the author of many publications, including several AIME 
technical papers. Fraser was born at Brockville, Ont., in 1885, 
and received the degree of BSc at Queen’s University. 


JAMES BAIN MORROW, vice-president of Pittsburgh Con 
solidated Coal Co., was awarded the Erskine Ramsey Gold 
Medal at the 1950 AIME Annual Meeting. He was born in 
Halifax, Nova Scotia, in 1886, and studied mining engineering 
at Dalhousie University. He has figured prominently on many 
committees of the AIME, and has been chairman, vice-chair- 
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man and Executive Committee member of the Coal Division. 
Many of his technical papers on the cleaning, sampling and 
preparation of coal have been published by AIME. In 1944 
Nichols Award. 


he received the Percy 


FOR DIRECTORS 


C. HARRY BENEDICT served as a director of AIME 
from 1945 to 1950 and as vice-president from 1948 to 1950. 
He was born at Pittsburgh, Pa., in 1876 and graduated from 
Cornell University in 1897. Immediately after graduation he 
became affiliated with Calumet and Hecla Consolidated Cop- 
per Co., and remained with this company until his retirement 
in 1948. Benedict developed the ammonia leaching process for 
recovery of copper from waste sand (tailing). He is writing 
a history of the company which will be published late this 
year. 

RUSSELL B. CAPLES was graduated from the Missouri 
School of Mines in 1910. Caples joined the Federal Min- 
ing and Smelting Co. in Idaho after graduation. He soon 
entered the employ of the Anaconda Copper Mining Co., where 
he assisted in the development of the electrolytic zinc process 
and the basic process for the Anaconda phosphate plant. In 
1918 he became superintendent of the Great Falls, Mont., 
plant., and was named manager in 194]. He has served on a 
number of committees of the Extractive Metallurgy Division 
of AIME, and is now a member of the Advisory Committee on 
Zine of the Defense Minerals Administration. 

THEODORE BENTON COUNSELMAN received the Engi- 
neer of Mines degree from Columbia University. He soon 
entered the ore dressing industry, and worked with various 
companies. Entering the equipment sales field, he joined The 
Dorr Co. as manager of the Central Division in Chicago. He 
is now in charge of exploiting the company’s newest develop- 
ment. Counselman has been chairman of the Chicago and 
New York sections, chairman of the Milling Methods Com- 
mittee, vice-chairman of the Blast Furnace Ceke Oven and 
Raw Materials Committee of the Lron and Steel Division, and 
has served on a number of other committees. He is currently 
chairman of the national Membership Committee. 

THOMAS L. JOSEPH is now chairman of the Iron and 
Steel Division. He born in Adamsville, Utah, in 1894, 
and received an MA degree from the University of Utah. He 
S. Bureau of Mines’ North Central 
Experiment Station in Minneapolis, and in 1936 became head 
of the Department of Metallurgy, Minnesota School of Mines 
and Metallurgy. where he is now assistant dean of the Insti- 
tute of Technology. During World War II he was a member 
of the War Metallurgy Committee. He received the J. E. John- 
son, Jr., Award in 1927 and the Richard W. Hunt Medal and 
Prize in 1937, and was Howe Lecturer for 1946. 

@ARL ERNEST REISTLE, JR., received a BS degree in 
chemical from the University of Oklahoma 
in 1922. He worked as a petroleum chemist for the U. S. 
Bureau of Mines at Bartlesville, Okla., and Laramie, Wyo. In 
1933 he became chairman of the East Texas Engineering 
Association, aiding in the establishment of conservation regu- 
lations. Since 1936 he has been with Humble Oil and Refining 
Co., and is now a member of the board and manager of the 
Production Department, located in Houston, Tex. During 
World War II he served as technical advisor for District ILI 
of the Petroleum Administration for War, and as national 
vice-chairman of the Oil Industry Advisory Committee for 
the OPA. x * * 


was 


was employed by the U. 


engineering 


SECTION 2 








Employment Notices 








The Jou RNAL will notices of 
men and jobs available. Companies and 
AIME members are invited to use this 
space, for no charge. 
Except as noted below, address replies 


post 


which there is 
to: Code (appropriate number), Jour- 
NAL OF PerroLteuM TecHNOLOGY, 408 
Trinity Universal Bldg., Dallas 1. Show 
return address on envelope. These re- 
plies will be forwarded unopened and 


no fees are involved. 


Replies to the positions coded Y5532, 
Y5531. and Y5480 below should be ad- 
Engineering Societies Per- 
sonnel Service, 8 West 40th St.. New 
York 18, N. , # The ESPs. on whose 
behalf these notices are published here, 
collects a fee from applicants actually 


dressed to: 


placed. 


POSITIONS 


@ Engineers (a) Petroleum engineer 
for field operations, 25-35, with degree 
in petroleum engineering or equivalent 
and about three to five years’ experi- 
ence. (c) Division safety supervisor, 
30-45, graduate, preferably mechanical, 
with five years’ experience in petroleum, 
super- 


drilling and production safety 


(d) Junior petro- 


petroleum 


vision. engineers 


leum, recent engineering 
graduates, capable of developing under 
training. Location, Peru, S. A. Y5532. 

@ Engineers. (a) Equipment and cost 
engineer, about 30, with petroleum en- 
gineering degree and about five years 
experience. (b) Production engineer 

5-40. with about four years” 


> y 2@ 95. 
Reservoir ¢ ngineer, 20 


petroleum 

experience. (c) 

10, with about five years” 

gineering experience. (e) Drilling engi- 
95.35 

a OOD, 


reservoir en- 


neer, with degree in petroleum 


engineering and about three years’ ex- 


perience. Location, Colombia, S. A. 


Y5531. 

@ Sales engineer, about 30, mechanical 
graduate, with some sales experience, 
for sale of industrial oil to paper mills, 
automotive manufacturers, heavy indus 
try, ete.. for major oil company. Loca 


tion, Wisconsin. Y5480. 


@ Applications now being received from 
petroleum, chemical or mechanical en 
gineers with PhD degrees and _petro- 


leum reservoir engineering experience 


for teaching positions on the petroleum 
engineering staff of a Southwestern in- 
stitution. Salary based on background. 


Initiative essential. Interviews during 


July and September. 
fective September A 


Appointments ef- 
Consulting and re- 
search encouraged. Send recent photo 


graph with data. Code 53 
PERSONNE! 


@ Engineer, petroleum production, reg- 
istered, presently employed, desires re- 
sponsible position with aggressive inde- 
pendent or small producing company 
Twelve years” experience mid-continent 
area all phases of production, special 
married 


izing water flooding. Age 32. 


Code 150. 


@ August graduate in petroleum engi- 
neering. MS degree. Age 22 
August 10 


two 


. single, non 
for 


-engi- 
non-engi 


available em- 


About 
neering field and office experience. Will 


veteran, 


ployment. years’ 


work anywhere. Code 151. 


@ Available 


duction superintendent; 20 


petroleum engineer - pro 


years’ €x- 
perience in drilling, production and sec- 


ondary recovery. Code 147. * * * 


Economic Aspects 

Continucd from Page 16, Section 1] 
force corporations to withhold 20 per 
cent of all wages, and pay that amount 
withheld direct to the government as 
taxes, 

One of the principal attacks in the 
general trend of nationalization or so- 
cialization has been against the oil and 
gas industry, for through the control 
of fuel and energy everything else can 
be controlled. Once private industry has 
lost control of oil and gas, the private 
control of other industries must like- 
wise eventually be lost. 

The position of oil in world affairs 
s not apt to be disturbed except by 
unwise legislation, but by such legisla- 
tion the capital of the oil industry can 
easily be taxed away and the present 
bountiful supply reduced to an actual 
shortage which thereby might well per- 
mit the government to take over. 

The chief weapon of attack on the 
oil business has been through the elim- 
ination or reduction of the depletion 
allowance granted to the oil industry 
before determination of the net income 
subject to taxes. In the present session 
of Congress, the Treasury Department 
sought to reduce the present rate of 
depletion, which is 274% per cent of the 
gross income, to 15 per cent, and previ- 
ously the Treasury had endeavored to 
have depletion eliminated entirely. For- 
tunately, due to the understanding of 
Congress of the subject, the Treasury 
Department has so far not prevailed. 
but the attack continues. 

Oil, when once found, produced and 
consumed, cannot be replaced. Basically, 
liquid energy and must be 


it is only 


in order to be 
~ «* * 


consumed or destroyed 


utilized. 





Proposed for Membership, Petroleum Branch 





Toto! AIME membership on January 31, 1951, 
was 17,228; in addition, 4,413 Student Associates 
were enrolled 

ADMISSIONS COMMITTEE 

Thomas G. Moore, Chairman; Carroll A. Gar 
ner, Vice-Chairman; George 8. Corless, F. W 
Hanson, Albert J. Phillips, Lloyd C. Gibson, R. D 
Mollison, John T. Sherman. Alternates, A 
Brinker, H. W. Hitzrot, Plato Malozemoff 
Given, T. D. Jones, W. H. Farrand 

Institute members are urged to review this list 
as soon as the issue is received and immediately 
to wire the Secretary's office, night message 
collect, if objection is offered to the admission of 
any applicant. Details of the objection should 
follow by air mail. The Institute desires to extend 
privileges but does not desire to admit persons 
unless they are qualified. Objections on appli 
cations should be received on the 15th of the 
month following publication in PETROLEUM TECH 
NOLOGY 

In the following list C/S$ 
status; R, reinstatement; M, 
Member; AM, Associate Member 
ciate; F, Junior Foreign Affiliate 


Ivan 


meons change of 
member; J, Junior 
S, Student Asso 


SECTION 2 


CALIFORNIA 
Banning Hansen, Beauford V 


Dunasky, John | J 
Barbarick, Ha 


ding (J 
Lee W M 


Riverdale 
COLORADO 
Loveland Weaver 
ILLINOIS 
Effingham 
. Charles 


KANSAS 
Wichita 
S-J); Sells, Harold R. (M 
LOUISIANA 

Harvey Trueblood, Harry 
Houma — Sylvester, Burke 
Lafayette Johnston, Richard E. (M 
Lake Charles Cummings, Aaron (J) 
New Iberia Kizer, Curtis A. (A 
OKLAHOMA 


Bartlesville 


Gibson, 


Lewis K 


Burke 


Sheldon 


Billy F J 
John A 


McGuire, Edward M. (J 
C/S-A-M 


Jack A. (M 


Dolman 
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PENNSYLVANIA 
ee Sittig, William R 
. Cresson 


Bartlebaugh, Dale L 


Jy 

State College Calhoun, John C. (M); 
Trevor W. (J). 
TEXAS 

Corpus Christi Whitmire, William A. 

Dallas Grubb, William E. (J). 

Fort Worth Acklin, Floyd, Jr. (J); 
lar, Wiiliam L. (J) (R, C/S-S-J) 

Houston Brooks, Aubrey J., Jr. (J) (C/S- 
S-J); Ray, Millard L. (M); Yust, Charles S 
(M; (C/S-A-M)} 

Levelland 
J-M) 
Lubbock 
McKinney 


(J). 
Tins- 


Edwards, Floyd C. (M) (C/S- 
Buckner, Roy E. (J) (C/S-S-J). 
Abbott, Seth J. (J) 
Losher, Don J. (J) (C/S-8-J) 
Stout, Kenneth R. (J). 
ARABIA 
Dhahran — Conroy, John Patrick (J) (C/S- 
S-J) **«rk 
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Petroleum Technology 
Adds Representative 


The firm of William R. Stewart, pub- 
lisher’s representatives, has been re- 
tained by the Petroleum Branch to act 
in the Midwest area in the solicitation 
of advertising for the JouRNAL oF Pr- 
rROLEUM TECHNOLOGY. The Stewart or- 
ganization will operate for PETROLEUM 
TecHNOLOGY in Illinois, Wisconsin, 
Michigan, Indiana, Ohio, Kentucky, 
Minnesota, lowa, Missouri, Kansas and 
Nebraska. Main offices of the Stewart 
Co. are at 9 S. Clinton St., Chicago, Il. 


With the conclusion of this transac- 
tion, representation for the JOURNAL OF 
PETROLEUM TECHNOLOGY is now com- 
plete in all major petroleum equipment 
manufacturing areas of the United 
States. The magazine is represented in 
the Southwest and West by the firm of 
McDonald-Thompson and in the East by 
Allen Associates. 


At the present, advertising on con- 
tract for 1951 amounts to a little more 
than $29,000.00. For the present ac- 
counting the publication is budgeted for 
$27.000.00 income from advertising. It 


is estimated that income from adver- 
tising will exceed $30,000.00 for 1951, 
an increase of about 60 per cent over 


similar income 1950. x * * 


Manpower Commission 
Issues First Bulletin 


Bulletin 1 of the Engineering Man 
power Commission of Engineers Joint 
Council, entitled “Utilizing Engineering 
Manpower Deferment Procedures,” 
was issued recently. It is a 14-page 
pamphlet, discussing the selection and 
utilization of technical personnel, spe 
cial Selective Service procedure per 
taining to college students and current 
college graduates, general occupational 
deferment procedures under Selective 
Service for non-reservists, recall of re- 
servists to active duty, and giving other 
pertinent information for engineering 
employers and employes. Copies have 
been supplied to all Local Section Sec 
retaries, or will be supplied at 25c each 
on individual order to Engineering Man- 
power Commission, llth Floor, 29 W. 
39th St.. New York 18,N. Y.*® *® *® 


AGI Reproduces Article 


On Conserving Manpower 


The American Geological Institute 
has made available in the public inter- 
est reprints of the article, “Manpower 
The Most Important Factor in Defense 
Mobilization,” which appeared in the 
May issue of the JournNAL or PeEtro- 
LEUM TecHNnoLocy. The article pre- 
sents a basic formula for conserving 
manpower in the best interests of both 
the military and industry. It reports on 
the current policy of the armed forces 
concerning selectees and reserve per- 
sonnel, and outlines a procedure for 
industry to follow in obtaining defer- 
ments for essential employees. The 
author, E. W. Berlin, is chairman of 
the API Defense Services Committee, a 
member of the National Petroleum 
Council's Manpower Committee, and a 
colonel in the U. S. Air Force Reserve. 
Copies of the article may be obtained 
from the American Geological Institute, 
2101 Constitution Avenue. Washington 
25, D. C. A limited supply is also avail- 
able in the office of the Petroleum 
Branch, AIME, Dallas. 7m ® 
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Book Reviews 











Marine Geology 

By Ph. H. Kuenen, John Wiley & a Inc., New York, 1950. 

551 pp., charts, tables, 914x6' in., cloth, $7.50. 
This book summarizes geological matter with regard to sub- 

marine landforms and their interpretation, sedimentary cover 

floor, the structure of atolls, etc. Present sedi- 

mentation processes and the evidence of conditions during 


diagrs., maps, 


of the sea 


accumulation that may be expected in ancient sediments are 
presented. Emphasis is placed on problems and relations rather 
than descriptions of properties, distributions, and numerical 
data. The author gives special attention to coral reefs, sub- 
marine canyons and the Moluccan deep-sea troughs, as topics 
with which he is especially familiar. 


Unit Operations 

By George Granger Brown and Associates, John Wiley & Sons, 
Inc., New York, 1950. 611 pp., illus., diagrs., charts, tables, 
1114x834 in., cloth, $7.50. 

This book is a comprehensive treatment of modern process 
operations, stressing the similarity of their basic principles, 
in order to clarify the problems of design and operation of 
equipment. A nomenclature is 
throughout in similar theoretical discussions. The operations 
are divided into four groups: fluids, separation by 
and rates of energy and mass 
transfer. Large illustrations and charts, as well as derivation 
of theoretical equations and descriptions of important equip- 
ment, the book’s usefulness. x~* 
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: LUFKIN FOUNDRY & MACHINE COMPANY 
LUFKIN, TEXAS 


Branch sales and service Houston, Dallas, New York, Tulsa, Los Angeles, Seminole, Oklahoma City 
-orpus Christ Odessa, Kilgore Wichita Falls) Casper, Wyoming; Great Bend, Kansas 





Lufkin Equipment in CANADA is handled by 
THE LUFKIN MACHINE co., LTD. 14321 108th Avenue, Edmonton, Alberta, Canada 
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